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Abstract 

While spike trains are obviously not band-limited, the theory of super-resolution 
tells us that perfect recovery of unknown spike locations and weights from low-pass 
Fourier transform measurements is possible provided that the minimum spacing. A, 
between spikes is not too small. Specifically, for a measurement cutoff frequency of 
fc, Donoho [2] showed that exact recovery is possible if the spikes (on M) lie on a 
lattice and A > 1/fc, but does not specify a corresponding recovery method. Candes 
and Fernandez-Granda [3,4] provide a convex programming method for the recovery 
of periodic spike trains (i.e., spike trains on the torus T), which succeeds provably 
if A > 2/fc and fc ^ 128 or if A > 1.26//c and fc ^ 10^, and does not need the 
spikes within the fundamental period to lie on a lattice. In this paper, we develop a 
theory of super-resolution from short-time Fourier transform (STFT) measurements. 
Specifically, we present a recovery method similar in spirit to the one in [3] for pure 
Fourier measurements. For a STFT Gaussian window function of width a = l/(4/c) 
this method succeeds provably if A > 1/fc, without restrictions on fc- Our theory 
is based on a measure-theoretic formulation of the recovery problem, which leads to 
considerable generality in the sense of the results being grid-free and applying to spike 
trains on both M and T. The case of spike trains on M comes with significant technical 
challenges. For recovery of spike trains on T we prove that the correct solution can be 
approximated—in weak-* topology—by solving a sequence of finite-dimensional convex 
programming problems. 

Keywords Super-resolution, sparsity, inverse problems in measure spaces, short-time Fourier 
transform 

Part of the material in this paper was presented at the IEEE International Conference on Acoustics, 
Speech, and Signal Processing (ICASSP), Florence, Italy, May 2014 [I]. 
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1 Introduction 

The recovery of spike trains with unknown spike locations and weights from low-pass Fourier 
measurements, commonly referred to as super-resolution, has been a topic of long-standing 
interest [5-12], with recent focus on f'^-minimization-based recovery techniques [3,13]. It 
was recognized in [2, 3,14-16] that a measure-theoretic formulation of the super-resolution 
problem in continuous time not only leads to a clean mathematical framework, but also 
to results that are “grid-free” [13], that is, the spike locations are not required to he on a 
lattice. This theory is inspired by Beurling’s seminal work on the “balayage of measures” in 
Fourier transforms [7,8] and on interpolation using entire functions of exponential type [9]. 
Specifically, de Castro and Gamboa [15] and Candes and Fernandez-Granda [3] propose to 
recover a periodic discrete measure (modeling the spike train), that is, a measure on the torus 
T, from low-pass Fourier measurements by solving a total variation minimization problem. 
Despite its infinite-dimensional nature this optimization problem can be solved explicitly, as 
described in [3,14,16]. Concretely, it is shown in [3,15] that the analysis of the Fenchel dual 
problem leads to an interpolation problem, which can be solved explicitly provided that the 
elements in the support set of the discrete measure to be recovered are separated by at least 
2//c, where fc is the cutoff frequency of the low-pass measurements, and fc ^ 128. More 
recently, Fernandez-Granda [4] improved the minimum distance condition to A > 1.26/fc, 
but had to impose the additional constraint fc ^ 10^. Donoho [2] considers the recovery 
of a spike train on M and proves that a separation of 1/fc is sufficient for perfect recovery 
provided that the spikes lie on a lattice; a concrete method for reconstructing the measure 
is, however, not provided. In a different context, Kahane [17] showed that recovery of spike 
trains on M under a lattice constraint can be accomplished by solving an infinite-dimensional 
minimal extrapolation problem, provided that the minimum separation between spikes is at 
least j-^y\og{l/fc). In [11,12] the recovery of periodic spike trains is considered in the context 
of sampling of signals with finite rate of innovation. The main result in [11] states that a 
periodic spike train with K spikes per period can be recovered from 2K + 1 Fourier series 
coefficients without imposing a minimum separation condition. The corresponding recovery 
algorithm falls into the category of subspace-based methods such as, e.g.. Multiple Signal 
Glassification (MUSIG) [18] and Estimation of Signal Parameters via Rotational Invariance 
Techniques (ESPRIT) [19], algorithms widely used for direction-of-arrival estimation in array 
processing. 
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1.1 Contributions 


In practical applications signals are often partitioned into short time segments and windowed 
for acquisition (as done, e.g., in spectrum analyzers) so that one has access to windowed 
Fourier transform, i.e., short-time Fourier transform (STFT), measurements only. Moreover, 
the frequency characteristics of the spike train to be recovered often vary over time, i.e., 
the spike locations can be more packed in certain time intervals, as illustrated in Fig. 1. 
Time-localized spectral information, as provided by the STFT, can therefore be expected 
to lead to improved reconstruction quality for the same measurement band-limitation. This 
motivates the development of a theory of super-resolution from STFT measurements, which 
is the goal of the present paper. Inspired by [2,3,14-16,20], we consider the continuous time 
case and employ a measure-theoretic formulation of the recovery problem. Our main result 
shows that exact recovery through total variation minimization is possible, for a Gaussian 
STFT window function of width a = l/(4/c), provided that the minimum spacing between 
spikes, i.e., the elements in the support set of the discrete measure to be recovered, exceeds 
l//c- While our theory applies to general window functions—from the Schwartz-Bruhat 
space—that are extendable to entire functions, the recovery condition A > 1 //^ is obtained 
by particularizing to a Gaussian window function of width a = l/(4/c). Similar recovery 
thresholds can be obtained for other window functions and for Gaussian window functions 
of different widths, but this would require adapting the computational parts of our proofs, 
in particular Appendices A and B. Our theory applies to spike trains on M and to periodic 
spike trains, i.e., spike trains on T, and we do not need to impose a lattice constraint on 
the spike locations. The case of general (i.e., without a lattice constraint) spike trains on 
M comes with signihcant technical challenges. For recovery of spike trains on T we prove 
that the correct solution can be approximated—in weak-* topology—by solving a sequence 
of hnite-dimensional convex programming problems. We hnally present numerical results, 
which demonstrate an improvement for recovery from STFT measurements over recovery 
from pure Fourier measurements in the sense of the minimum spacing between spikes being 
allowed to be smaller for the same cutoff frequency. This improvement comes, however, at 
the cost of increased computational complexity owing to the redundancy in the STFT, which 
leads to an increased number of measurements and thereby a larger optimization problem 
size. 

1.2 Notation and preparatory material 

The complex conjugate of G C is denoted by The hrst and second derivatives of the 
function ip are designated by p' and p", respectively, for n G hl\{ 0 , 1 , 2 }, its nth derivative 
is written as p^'^\ The sine function is defined as sinc(t) := sin(t)/f, t 7 ^ 0, and sinc(O) := 1. 
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Figure 1: Example of a spike train. The spikes are more densely packed in certain intervals 
than in others. 


(5o/ is the Kronecker delta, i.e., = 1 for £ = 0 and = 0, else. Uppercase boldface 

letters stand for matrices. The entry in the kth row and Uh column of the matrix M is 
mfc/. The superscript ^ denotes Hermitian transposition. We dehne the real inner product 
of the matrices X, Y G as (X, Y) := Re{Tr(Y^X)}. For a hnite or countable set 

denotes the space of bounded sequences a = with corresponding norm 

llall^oo = sup^gf^ \a(\. Linear operators are designated by uppercase calligraphic letters. For 
n G M, 7^ denotes the translation operator, i.e., (Tu^)(t) := ip(t — u), for all t G M. For 
/ G M, Wl/ is the modulation operator, i.e., {Aifip)(t) := for all t G M. Let X 

and Y be topological vector spaces, and X* and Y* their topological duals. The adjoint 
of the linear operator £; X —)■ Y is denoted by C* \ Y* —)■ X*. For an extended-valued 
function <^9 : X —)■ M U {cxd}, we use dom(p to denote its domain, i.e., the subset of X where 
(p takes hnite value, and cont p stands for the subset of X where p takes hnite value and is 
continuous. The function p: X —)■ M U {c)o} is lower semi-continuous if, for every a G M, 
the set {x G X: p{x) ^ a} is closed. Let {x,x*) be a dual pairing between x G X and 
X* G X*. \i p\ X —)• M U {c)o} is a convex function, its Fenchel convex conjugate is the 
function (p*: X* —)■ M U {oo} dehned by p*{x*) := sup 3 ,g^ {(x, x*) — p{x)}^ for all x* G X*. 
If p is not identically equal to oo and if xq G X, dpi^xo) denotes the subdiherential of p 
at the point xq, i.e., dp^Xo) := {x* G X*: p{x) ^ p{xo) + {x — Xo,x*), for all x G X}. The 
set of all solutions of an optimization problem (P) is denoted by Sol{(P)}. For a measure 
space (X, and a measurable function (p: X —)■ C, the integral of p with respect to the 
measure /i is written as p{x)p{x), where we set dx := A(x) if A is the Lebesgue measure. 
For p G [1, cxd), L^(X, E,p.) denotes the space of measurable functions <p: X —)■ C such that 
II^IIlp := (Jx < OO- The space L°°(X, S,/i) contains all measurable functions 

<p: X —)■ C such that ||</9 ||l°° •= inf{C' > 0: |</9(x)| ^ C, for //-almost all x G X} < oo. 
For p G LP(X,E, fi) and ip G L'^(X,S,/i) with p,q G [1, C)o] satisfying 1/p-|- 1/q = 1, 
we dehne the complex inner product (plpj) ■= p(t)'ip(t)p{t) and the real inner product 
{p, p)) := Re I (p(x)'?/(x)/i(x)|. For a separable locally compact Abelian topological group 
G (e.g., the additive group M or the torus T := M/Z, both endowed with the natural topology), 
we write L'^{G) in the particular case where S = B{G) is the Borel cx-algebra of G and p 
the Haar measure on G. S{G) is the space of Schwartz-Bruhat functions [21]. We will not 
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need the entire formalism of Schwartz-Bruhat functions, as we exclusively consider the cases 
(7 = 1^ and G = T. Specihcally, iS(M) is the space of Schwartz functions, i.e., functions 
9 ?: M —)■ C that are inhnitely often differentiable and satisfy supjg]^ |t™| < oo, for all 

n,m E N. iS(T) is the space of inhnitely often differentiable functions. A complex-valued 
bounded hnitely additive measure p on G is a function p: i3(G) —)■ C such that for all disjoint 
hnite collections of sets in B{G), 





( 1 ) 


e=i 


and for every B G B{G) we have \fi{B)\ < oo. A complex-valued bounded hnitely additive 
measure is said to be regular if for each B G B{G) and each e > 0, there exists a closed set 
FEB and an open set G ^ i? such that for every G E G \ F, |/i(G)| < e. We denote 
the space of comp lex-valued regular bounded hnitely additive measures on G by Ai{G). For 
t E G, St E Ai{G) designates the Dirac measure at t, dehned as follows: for B E B{G), 
St{B) = 1, if t G B, and St{B) = 0, else. The support supp(p) of p G Ai{G) is the largest 
closed set G C G such that for every open set B E B{G) satisfying i? fl G 7 ^ 0, it holds that 
jj,{B nG) 7^ 0. We dehne the total variation (TV) of /i G Ai{G) as the measure |/i| satisfying 

VBee(G), 71(B) := sup W 7(^)1, 

where 11(5) denotes the set of all hnite disjoint partitions of 5, i.e., the set of all disjoint hnite 
collections of sets in B(K) such that B = ljt=i Throughout the paper, we equip 

the space Ai{G) with the TV norm ||/7||ty := \n\{G). By [22, Thm. IV.6.2], M(G) is the dual 
of the space Gf,(G) of complex-valued bounded continuous functions (f. Gb(G) is equipped 
with the supremum norm = sup^g^ |(p(t)|. By analogy with the real inner product 

in T^(G), we dehne the complex and real dual pairing of the measure ju E M(G) and the 
function (f E Gb(G) as (^^Ip) := and ( 93 , p) := Re | respectively. We 

endow M(G) with the weak-* topology [23, Chap. 3], i.e., the coarsest topology on Ai{G) 
for which every linear functional Ai{G) —)■ M dehned by p h-)■ = {‘-P-ijj), with 

(p E Gf,(G), is continuous. The Pontryagin dual group of G is denoted as G, and the Fourier 

transform of /i G Ai{G) is the function p: G —)■ C dehned by /!(/) = / for all 

f E G. If G is equipped with the metric I'l, we denote by Br{G) the ball in G — with respect 
to the metric I'l—that is centered at 0 and has radius r, i.e., Br{G) := {/ G G: |/| ^ r}. 
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2 The problem statement 


We first state the recovery problem considered and then discuss its relation to prior work. 
Let G = M or G = T (for G = M, we have G = M, and for G = T, we get G = Z). We model 
the spike train with weight G C\{0} attached to the point G G G by a measure in M.{G) 
of the form 

^ = ( 2 ) 

where hi is a hnite or countably inhnite index set. The measure is supported on the set 
T := assumed closed and uniformly discrete, i.e., there exists a <5 > 0 such that 

\t(. — tf'l ^ 5, for all G n. We have ||/i||Ty = I®^l- Moreover, since fi G Ai{G), we 

also have ||/i||rpy < oo. Henceforth /i exclusively designates the measure dehned in (2). 
Suppose we have measurements of jj, in the time-frequency domain of the form 

i/(g /) = (Vg/i)(r, /), T eG, f e Bf^{G), 

where fc is the cutoff frequency (when G = T, the cutoff frequency fc becomes an integer, 
which we denote by K) and 

(Vg/i)(r,/) := [ c/(t - (3) 

Jg 

denotes the STFT [24] of jj, with respect to the window function g taken to be a Schwartz- 
Bruhat function which, in addition, is assumed to be extendable to an entire function (exam¬ 
ples of such functions are the Gaussian function, Hermite functions, and the Fourier transform 
of smooth functions with compact support). 

We are interested in recovering the unknown measure g, from the STFT measurements y 
through the following optimization problem: 


(STFT-SR) minirnize ||zz||rpy subject to y = AgU, 
where Ag: M.{G) —)■ L^(G x G) maps v G M.{G) to the function p G L^(G x G) given by 


V(r,/)GGxG, p(r,/) 


(V,z/)(r,/), /G%(G) 

0 , otherwise. 


(4) 


The idea of minimizing the TV norm to recover g from y originates from the seminal work 
of Beurling [9], who studied so-called “minimal extrapolation”. The reason for hoping that 
(STFT-SR) delivers the correct solution lies in an observation made in [3,9,13,14,16], which 
states that the TV norm Il'Urpy acts as an atomic-norm regularizer as a consequence of the 
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extreme points of the unit ball {u G Ai{G): ||z^||ty ^ 1} being given by the Diracs Sx, 
X & G. Minimizing H-Ilrpy therefore forces the solution to be discrete much in the same 
way as minimizing the f'^-norm for vectors in enforces sparsity. For more details we 
refer to [14], which contains a general and rigorous analysis of inverse problems in spaces of 
measures. 


3 Relation to previous work 

Super-resolution theory dates back to the pioneering work by Logan [5,6] and by Donoho [2]. 
Specifically, [ 2 ] considers recovery of the complex measure 

h = X] (5) 

t&n 

in A1(M), where ae G C\{0} and tg G aZ, for all i E fl, with a > 0, from pure Fourier 
measurements 

y(f) = A(/) := / / € [-/../J, 

Jr 

where fc is the measurement cutoff frequency. The main result in [2] is as follows. If the 
support supp(/r) = of /i has Beurling density 

r-.+/ / I- n+(supp(/i),r) . 

Il+(supp(/i)) := hmsup- < fc, 

r^oo 

where for r > 0, n+(supp(/i), r) denotes the largest number of points of supp(/i) contained in 
the translates of [0,r], then y uniquely characterizes y among all discrete complex measures 
z/ G A1(M) of support supp(z/) C aZ with Beurling density strictly less than fc. 

In [3] Candes and Fernandez-Granda deal with the recovery of periodic spike trains of 
the form (5). The resulting measure y = ^1,^=1 with ti G T, for all i E {1, 2,... , L}, 
is in A^(T) (i.e., G = T), and the goal is to recover y from the band-limited pure Fourier 
measurements 

yk = fi{k), k E (6) 

where K corresponds to the (integer) cutoff frequency. This is accomplished by solving the 
following optimization problem: 

(SR) minimize ||zz|L,, subject to y = JFv, 


where y := {yk}k=-K ^ and JF: A^(T) —)■ maps z/ G A^(T) to the vector 



X {xk}k=-K e (j^ 2 _ft:+i 2 ;^ := z/(/c). The main result in [3] establishes that /x is the 

unique solution of the problem (SR) if the minimum wrap-around distance 

A := min min + ^| 

between points in T = obeys A > ‘IjK. Moreover, it is shown in [3] that, under 

A > 2/7^, the support T = of p can be recovered by solving the dual problem 

(D-SR) maximize (c, y) subject to ||J^*c|| ^ 1 

ceC^-ff+i 

of (SR), where the adjoint operator J^* is given by 

K 

Vc e 

k=-K 

The dual problem (D-SR) is equivalent to a hnite-dimensional semi-dehnite program. The 
tg are determined by locating the roots of a polynomial of hnite degree on the unit circle, 
built from a solution of (D-SR). More recently, Fernandez-Granda [4] improved the minimum 
distance condition to A > 1.26//c, but had to impose the additional constraint fc ^ 10^. 

In summary, while Donoho considers super-resolution of measures on M with possibly 
inhnitely many ti to be recovered and hnds that A > 1//^ is sufficient for exact recovery, he 
imposes a lattice constraint on the locations of the t£ and does not provide a concrete recovery 
method. Candes and Fernandez-Granda [3,4] on the other hand provide a concrete recovery 
method which succeeds provably for A > 2//c and fc ^ 128 or A > 1.26/fc and fc ^ 10^, 
but is formulated for measures on T only and hence needs the number of unknown locations 
tg to be hnite. The main contribution of the present paper is to show that recovery from 
STFT measurements through convex programming succeeds for A > l//c, both for measures 
on M and T, without imposing a lattice constraint and without additional assumptions on 
fc- The rigorous treatment of the case of measures on M comes with signihcant technical 
challenges as the number of atoms of the measure to be recovered can be inhnite and no 
lattice constraint is imposed. 

A polynomial root-hnding algorithm for the recovery of the spike locations tg & T was 
presented in the context of sampling of signals with hnite rate of innovation [11]. Specihcally, 
the algorithm proposed in [11] can be cast as a subspace algorithm and does not need a min¬ 
imum spacing constraint; corresponding performance results for the noisy case were reported 
in [25,26]. 

Very recently [27] considered a super-resolution problem similar to (SR), aimed at the 
recovery of a complex Radon measure /i on T'^, d G N\{0}, from its Fourier coefficients /i(/c). 
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for /c G A, where A is an arbitrary finite subset of l/. The main result of [27] is an adaptation 
(to the torus) and a generalization to higher dimensions of Beurling’s theorem on minimal 
extrapolation. On a conceptual level [27] establishes a connection between Beurling’s theory 
of minimal extrapolation [9, Thm. 2] and the work by Candes and Fernandez-Granda on 
super-resolution [3,4], 

4 Reconstruction from full STFT measurements 

Before considering the problem of recovering /i from /c-band-limited measurements y, we 
need to convince ourselves that reconstruction is possible from STFT measurements with 
full frequency information, i.e., for /c = oo. 

It is well known that the STFT of a function ip G L^(]R) with respect to the (nonzero) 
window function g, defined as 

V(r,/) G {Vg(p){Tj) := {MfTrg\(p) = [ <^(t)c/(t - (7) 

Jr 

uniquely determines (p in the sense of Vg(p = 0 implying ip = Q. The signal p G L^(M) can be 
reconstructed from VgP via the following inversion formula [24]: 

= [0^g^){r,f)Mfrrg drdf. (8) 

\\g\\L2 Jr Jr 

Since for G = T, G = Z, the STFT of a function 99 G T^(T) is given by 

V(r, k) eT xZ, {Vgp){T, k) := {J^kTrg\p) = f p{t)g{t - 

Jo 

with the corresponding reconstruction formula 

1 V 

^ / Clos'd )(a k)MkTrg dr. 

Grochenig [24] extended the classical definition of the STFT for functions in L‘^{G) to 
tempered distributions. As measures in Ad(G) are special cases of tempered distributions, 
this extended dehnition applies to the present case as well. The remainder of this section 
is devoted to particularizing the STFT inversion formula [24, Gor. 11.2.7] to measures in 
A1(G). In the process, we obtain a simple and explicit inversion formula for discrete measures. 
Proposition 2, whose proof will be presented in detail as it provides intuition for the proof of 
our main result in the case of incomplete measurements. 

The STFT of a measure u G A1(G) is obtained by replacing the complex inner product 
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on L‘^{G) in (7) by the complex dual pairing between measures u G M(G) and functions 
(f E Gb(G), that is, '■= 

Definition 1 (STFT of measures). The STFT of the measure v E M.{G) with window 
function g E 5(G) \{0} is defined as 


(Vgz/)(r,/) :={MfTrgW)= [ 9{t-T)e 

Jg 

The following result provides a formula for the reconstruction of u E Ai(G) from Vgio. 


Proposition 1 (STFT inversion formula for measures, [24, Cor. 11.2.7]). Let v E M(G) and 
g E 5(G)\{0} he such that \\g\\i 2 = 1- The quantities 



(yg^) {'r,f)M fTrg drd/, for G = R, 


Y. I (VgZ/)(r, k)MkTrg dr, for G = T, 


fcez ' 


define measures in A1(M) and A1(T), respectively, in the sense that 


V<^ e a(G), 




[ [ (VgZ/)(r,/)((^|Af/7;^)drd/, 

M. M ^ 

f y9’^)i^J)iT\MkTrg)dT, 

fcez 


Moreover, 


V = 


[ /(VgZ^)(r,/)Af/7;^ drd/, 

»/ M »/ M ^ 

X] / yg^){^,k)MkTrg dr, 
fcez “^0 


G = M 
G = T. 


G = M 
G = T. 


(9) 


An immediate consequence of Proposition 1 is the injectivity of the STFT of measures. 
To see this, simply note that Vgio = 0 in (9) implies y\i') = 0 for all (p E Gb(G), which means 
that necessarily z/ = 0. Equivalently, for z/i, z/2 G M(G), as M.{G) is a vector space and hence 
z/i — z /2 G M.{G), it follows that Vg^i = VgV 2 with g E 5(G)\{0} necessarily implies z/i = z/ 2 . 
We can therefore conclude that every measure in Ai{G) is uniquely determined by its STFT. 
While the inversion formula in Proposition 1 applies to general measures in M(G), we can 
get a simplihed inversion formula for discrete measures, as considered here. Specihcally, for 
discrete measures inversion reduces to determining the support set T = and the 

corresponding complex weights only. 


Proposition 2 (STFT inversion formula for discrete measures). Let g E 5(G) \{0} be such 
that || 5'||^2 = 1- Let pi, ;= E M.{G), where is a finite or countably infinite index 
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set and ag G C\{0}, for all i E il. The measure /i can be recovered from its STFT Vgfv 
according to 


^ - r)e 2 "^Mrd/ 


l-F 

K „i 


lim ——- 

K^oo 2K + 1 


^ [ i^ 9 T){r,k)g{t- 

k=-K-^o 





t = tg 
otherwise 

t = tg 
otherwise 


G = R (10) 
G = T. (11) 


Proof. We prove the proposition for G = R and note that the proof for G = T is similar. Let 
t eR. For every F > 0, we have 

^ r /(V,fi)(T,/)9(«-r)e“^‘dTd/ 


l-F . 


[ [ - r)git - ‘"Mrd/ (12) 

J-F Jr 

9{te-r)g{t-T)dT] (^[ e‘^^^kit-te)^j\ ^^ 3 ^ 

g^n VJr / \ J-F J 

'^2agR{t — tg) smc{2FF{t — tg)), (14) 


2F 


e&n 


where R E iS(M) is the autocorrelation function of g, that is, 


Vt G M, R{t) := / g{T)g{t + r)dr. 


(15) 


Thanks to Fubini’s Theorem, the order of the integral on [—F, F] x M and the sum in (12) 
can be interchanged as, using the Cauchy-Schwarz inequality, we have 


ag 


^2F 


eefi 


l-F ■ 


g{tg - T)g{t - T)e^^^k{t h) ^ ||^|| 2 ^ ^ ^ 




since g E 5(M) C L‘^(R), and, by assumption, '^g^Q \ag\ < oo. Now we want to let F —)■ cx) 
in (14). First note that 


lim sinc(27rF(t — tg)) 

F^oo 


1 , t=tg 
0 , otherwise. 


For the limit F —)■ cx) of the series in (14) to equal the sum of the limits of the individual 
terms, we need to show that the sum converges uniformly in F. This can be accomplished 
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5 Reconstruction from incomplete measurements 

We are now ready to consider the main problem statement of this paper, namely the re¬ 
construction of fi from band-limited (i.e., fc < oo) STFT measurements via (STFT-SR). 
Henceforth, we assume fc < oo. 

We hrst verify that the range of the operator Ag dehned in (4) is indeed L^{G x G), and 
then show that Ag is bounded. 

Lemma 1 (Range of Ag). Let g G iS(G)\{0}. AgV G L^(G x G) for v G M.{fG). 

Proof. We provide the proof for G = M only, the proof for G = T again being essentially 
identical. Let u G W1(M). We have the following: 

[ [ \iA^)ir,f)\dTdf= j [ |(VgZ/)(r,/)|drd/ 

= / 1 / f clrt d/ 

J-fc \Js. / 

I If ( f A - l>'l(*)) dr] d/ 

= ‘Jfc [ ([ \g{t-T)\\u\{t)]dT. (16) 

JR \JR / 
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Since g G L^(M.), we have 


[ ( [ \9{t-T)\dT\\u\{t) = \\g\\^, ||z/||tv<oo. 

Jr \Jr / 

Fubini’s Theorem then implies that the integral in (16) is hnite, thereby concluding the 
proof. □ 

Lemma 2 . Let g G 15 ( 6 *)\{ 0 }. The operator Ag: M.{G) — )■ x G) is bounded. 

Proof. Again, we provide the proof for G = M only, the arguments for G = T are essentially 
identical. For every u G Ai{G), we have 




-fc 


^ f P P f 

IIA^IIli=/ / l(V<;Z^)(r,/)|drd/= 

~fc 

^ I / l^(f-'r)| |z/Kt)drd/= 2 /, 

'-fc Jtr 



= 2fc 



g{t-T)e 

\g{t-T)\ \ir\{t)dT 
|^(t-r)|dr|z/|(t) = 2/c||c/||^i ||z/||tv> 



(17) 


□ 


where we used Fubini’s Theorem in the hrst equality of (17). 

The adjoint of the operator Ag plays a crucial role in the ensuing developments. 

Lemma 3 (Adjoint of Ag). Let g G iS(G)\{0}. The adjoint operator A*g'. L°°{G x G) —)■ 
Gf,(G) of Ag maps c G L°°{G x G) to the funetion 


c(r, f)g{t - rje^^^-^Mrd/, 


’-fc JR 

fc «1 


t I — t ^ / c(r, f)g{t - 


G = M 
G = T. 


( 18 ) 


t—f. 


Moreover, it holds that A*g* = Ag. 

Proof. Again, we provide the proof for the case G = M only, the arguments for G = T being 
essentially identical. 

First, we note that L°°{Mf) is the topological dual of L^{Mf). Therefore, A*g maps L°°{Mf) 
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to (A^(M))*. For V G A^(M) and c G we have 

rF r _ 


{AgU, c) = Re 
= Re 

= Re 

= Re 


'-fc 

rfc 

'-fc 


{Agiy){Tj)c{Tj)dTdf 


c{r, f)g{t - r)e drd/ 


c(r, f)g{t - r)e ^"^Mrd/ u{t) 


'-fc 


(19) 


i: 


c{r, f)g{t - r)e2’^RMrd/j z/(t) 


= {iy,A*gc), 

where (19) follows from Fubini’s Theorem whose conditions are satisfied since 



|c(r, /)| \g{t - r)| drd/ \u\{t) ^ 2/^ ||c|| 


-fc 


Li Ir IItv ’ 


( 20 ) 


as a consequence of c G L°°(M^), G iS(M), and z/ G A1(M). 

It remains to show that A*gC G Cb{^) for all c G L°°(M^). To this end, we first note that 
for c G we have 


''fc 


VfGM, {A*gC){t)= / c(r,/)^(t - r)e^^*-^*drd/ 


'-/c 

'-fc 


c{t - u, f)g{u)e^^"^Audf. 


As g E 5(M), the mapping t ^ c(t — u, f)g{u)e^'^^A jg continuous on M for all {u, f) ^ 1^^- 
Moreover, we have 


VtGM, V(M,/)GM^ |c(t-M,/)^(M)e^^*'^*| ^ ||c||^^ |^(m)| 


and (m,/) e-)■ ||c||^oo \g{u)\ is Lebesgue-integrable on M x [—fcfd- R therefore follows, by 
application of the Lebesgue Dominated Convergence Theorem, that A*gC is continuous on M. 
Furthermore, we have 

VtGM, \{A*gC){t)\ ^ [ /" |c(r,/)| |^(t-r)|drd/ ^ 2/c||c||^„o ll^llii < cx), 

J—fc Jk. 

as c G L°°(M^) and g G iS(M). This shows that A*gC is bounded and therefore A*gC G Cf,(M), 
for all c G ^“(M^). 

We next show that A*g* = Ag. As M*: L°°(M^) —)■ Cf,(M) and A1(M) is the topological 


15 



dual of we have A*g* A^(M) —)■ Let v G A^(M) and note that 





= {c,A-;v) 


c(r, f)g{t - 


LMrd/^ 


-fcJ^ / 

(WT){i'At-T)e-^-''‘“v{t)] 

R \JR / 

At, f )(Agv)(T, /)dTd/| 



(21) 


where (21) follows by Fubini’s Theorem whose conditions are satished thanks to (20). This 
establishes that A**^ = AgV, for all z/ G A^(M). Since AgV G L^(M^), for all z/ G A^(M), it 
follows that A*g* A^(M) —)■ L^(M^), which completes the proof. □ 

Since the space M(G) is inhnite-dimensional, proving the existence of a solution of 
(STFT-SR) is a delicate problem. It turns out, however, that relying on the convexity 
of the TV norm ||■||^py and on the compactness—with respect to the weak-* topology—of the 
unit ball {z/ G M(G): ||z/||rpy ^ 1}, the following theorem ensures the existence of a solution 
of (STFT-SR). 


Theorem 4 (Existence of a solution of (STFT-SR)). Let p := aAte ^ Ai{G), where hi 
is a finite or countably infinite index set and G C\{0}, for all £ G Ll. Let g G iS(G)\{0}. 
If y = Agfi, there exists at least one solution of the problem (STFT-SR). 

Proof. As y = AgP the set iV := {z/ G Ai{G) : y = AgP} is non-empty, and hence. 


m := inf {||z/||rpy : z/ G N} < oo. 


For all n G hl\{0}, we dehne the intersection of N and the closed ball of radius m + 1/n in 
M(G) according to 

r „ „ 1 

Hn := < z/ G iV: z/ Ly ^ m H— 

[ n 

Since N is the preimage of {y} under the linear operator Ag, N is convex. Fnrther, as 
Ag is linear and bonnded, Ag is continnous with respect to the topology induced by the 
norm IHIrpy, which implies that N is closed with respect to the norm IHIxv Therefore, for 
every n E N \{0}, Kn is convex, closed with respect to the norm IHIrpy, and bonnded. By 
application of [23, Cor. 3.22], every Kn, n G hl\{0}, is compact with respect to the weak-* 
topology. Since the sets Kn, n G N\{0}, are non-empty and Kn' C Kn for all n,n' G N\{0} 
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such that n ^ n', by dehnition, we get 

OO 

Pi 7^ 0 . 

n=l 

Hence, we can hnd at least one z/q in which, by dehnition, satishes ||z^o|Itv ^ m + l/n 

for all n G hl\{0}. Letting n —)■ cx), we obtain ||z^o|Itv ^ Since m also satishes m ^ 
||z/o||ty, it follows that ||z^o|Itv ~ This shows that z/q is a solution of (STFT-SR), thereby 
completing the proof. □ 

Theorem 4 shows that (STFT-SR) has at least one solution. Next, with the help of 

Fenchel duality theory [28, Chap. 4], we derive necessary and sufficient conditions for /i to be 

contained in the set of solutions of (STFT-SR). Specihcally, we apply the following theorem. 

Theorem 5 (Fenchel duality, [28, Thms. 4.4.2 and 4.4.3]). Let X and Y be Banach spaces, 
99 : X — )■ M U {c)o} and R —)■ M U {c)o} convex functions, and C: X ^ Y a bounded linear 
operator. Then, the primal and dual values 

p = inf {(p{x) + '4){Cx)} (22) 

x£X 

d= sup {-ip*{C*y*) (23) 

y*eY* 

satisfy the weak duality inequality p ^ d. Furthermore, if ip, fj, and C satisfy the condition 


{C dom p) n cont t/’ 7^ 0, 


(24) 


then strong duality holds, i.e., p = d, and the supremum in the dual problem (23) is attained 
if d is finite ^. 


Application of Theorem 5 yields the following result. 

Theorem 6 (Fenchel predual). Let y := aeSt^ G M.{G), where hi is a finite or countably 
infinite index set and ai G C\{0}, for all i E fl. Let g G 5(G) \{0} and y = Agp. The 
Fenchel predual problem of (STFT-SR) is 


(PD-STFT-SR) maximize {c,y) subject to ||.4,*c|| 

ceL°-{Gxd) °° 

with the adjoint operator A* as in (18). In addition, the following equality holds 


mm < z/ 


TV : AgU = y,u e 

= sup{ {c, y) : I|.4jc||_^ < l,c e L“(G x G) 


(25) 


^Here, p,d gRU {±00}. 
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Moreover, if (PD-STFT-SR) has a solution cq G L°°{G x G), then 

IJ supp(z/o) C {t e G: |(^*Co)(t)| = 1}. (26) 

!^oeSol{(STFT-SR)} 

Proof. Similarly to what was done in [14, Prop. 2], we prove that (STFT-SR) is the dual of 
(PD-STFT-SR) by applying Theorem 5 with X = L°°{G xG),Y = Gb{G), C = A*g and with 
the functions 

VceL“(GxG), cpic):=-{c,y), (27) 

and 

Vr/ e Gfe(G), := xsiv) •= i ^ ^ ^ (28) 

I oo, otherwise, 

where Xb is the indicator function of the closed unit ball B := {rj E Gb{G): ||?7||oo ^ 1}. 
Thanks to Lemma 3, we have C* = A*g* = Ag. As (f is linear it is convex on L°°(G x G). The 
function ip is convex on Gb{G), because the closed unit ball B is convex and the indicator 
function of a convex set is convex. The Fenchel convex conjugate of ip is 

Ve e L\G X G), = sup |(c,0 - ^{c ): c G L“(G x G)} 

= sup I (c, -I- y) : c G (G X G) I 

= Xm(-0, (29) 

where the last equality follows from the fact that the set |(c, .^-|-|/) : c G L°°{G x G)| is 
not bounded when ^ —y. As for the convex conjugate of ip, it is known that the convex 
conjugate of the indicator function of B is the support function of B, that is, 

Vz/ G M(G), ip*(n) = CBin) := sup {(y, n) : y E Bj 

= sup {(y, n) :yE Gb{G), ||r 7 ||^ ^ 1} 

~ II^IItv ’ (80) 

where (30) follows from [22, Thm. IV.6.2]. It then follows by application of Theorem 5 that 
the dual problem to (PD-STFT-SR) is 

(STFT-SR) minimize ||zz||rpy subject to y = AgU. 

iy€M(G) 

Next, we prove (25). To this end, we verify (24), which implies (25) by strong duality, i.e., 
p = d, as follows. We start by noting that with X = L°°{G x G),Y = Gf,(G), C = A*g, and 
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9 ?, tjj as defined in (27), (28), we get from (22) that 


P= inf ^ i <^(c) + '4^{,A*c) \= inf ^ - (c, y) + Xb{A*c) 

c&L°°(GxG) ■’ c&L°^{GxG) 

= sup |(c,|/) - XB(^gC)| = sup|(c,|/) ; ||^*cL ^ l,c e L“(G'X G)| , 

ceL°°{GxG) 

where in the last step we used the fact that xb{A*c) is hnite, actually equal to zero, only for 
c G L°°{G X G) such that ||x4.*c||^ ^ 1. With (p*, ip* as in (29), (30), and (Cf,(G))* = A4(G), 
it follows from (23) that 

d= sup I — p*{Agiy) — ip*{—iy)\ = sup I — p*{—Agiy) — ip*{iy] 

,,c\A(n\ t J ueMiG) f 


lyeMiG) 

= sup 

ueM{G) 




ZZ 


ly 


TV 


= inf 

ueM(G) 


,v) + \\v\ 


TV 


= inf I llnq^, : A^v = y,v £ 

= min I ||n|L„ : A.v = y,i^€ 


(31) 

(32) 

Again, in (31) we used the fact that X{y}{Agiy) is hnite, actually equal to zero, only for 
ly G Ai{G) such that Agiy = y. The inhmum in (31) becomes a minimum in (32) as a 
consequence of the existence of a solution of (STFT-SR), as shown in Theorem 4; this solution 
will be denoted by pq below. It remains to establish (24). By dehnition of ip, we have 
dom(/i = L°°{G X G), and 0 G £dom((^). Furthermore, since cont'^ = R, we have 0 G cont'^ 
and therefore 0 G (/idomi/i) fl (cont'^), which shows that (24) is satished. 

Next, we verify (26). If (PD-STFT-SR) has at least one solution, say cq G L°°{G x G), 
then necessarily ||zzo||rpy = {co,y) as a consequence of strong duality. Since z/q is a solution of 
(STFT-SR), we have y = AgVo, and consequently. 


TV 


= (co,l/) = (co,x4gZ/o) = {A*gCQ, z/q) = Re M (^;co)(t)z/o(t) 


'G 


Applying [29, Thm. 6.12] to z/q, we can conclude that there exists a measurable function 
rp. [0, 27r) —)■ M such that 




TV 


= Re 


'G 


{A*co){t)iyo{t) [“ = Re <; / (^*co)(t)e*’'^*Vo|(f) = Re <( / h{t)\po\{t) )■, 


'G 


IG 
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where we set h := (^*co)e Writing h = |/i|e*®, where 9: [0, 27r) —)■ M, we then have 

0 = I^oIItv - Wo\{t) - Re |/i(t)|e"*®W|z/o|(^)| 

= [ (l - I^W|cos(0(t)))|z/o|(t). (33) 

V G- 

As L°°{G X G) —)■ Gfe(G), the function \h\ = |x4*co| is continuous. Moreover, since cq is 
a solution of (PD-STFT-SR), it follows that \\h\\^ = ||x4.*Co||^ ^ 1, and hence, 

yt e G, 1 — \h{t) \ cos( 6 '(t)) ^ 0. 

Therefore, (33) implies that \h(t) \ cos{9(t)) = 1 for |z/o|-almost all t E G, which, as a conse¬ 
quence of cos ^ 1 and ||/i||^ ^ 1, yields \h(t)\ = 1, for |z/o|-almost all t E G. We complete 
the proof by establishing that for all t E supp(z/o), \h{t)\ = |(M*co)(t)| = 1. Since \h{t)\ = 1 
for |z/o|-almost all t E G, the set S' = {t G G: \h(t)\ ^ 1} satishes |z^o|(>S') = 0. We need 
to show that S fl supp(z/o) = 0- To this end, suppose, by way of contradiction, that there 
exists a to E supp(z/o) such that \h(to)\ ^ 1. Then, S fl supp(z/o) 7 ^ 0 as to ^ “S' fl supp(z/o). 
Since \h\ is continuous on G and M\{1} is open in M, the set S is open. Hence, by dehnition 
of supp(z/o), we have |r^o|(S') ^ l^^oKS* fl supp(z/o)) > 0, which contradicts Iz^oKS*) = 0 and 
thereby completes the proof. □ 

We emphasize that (PD-STFT-SR) is the predual problem of (STFT-SR), i.e., (STFT-SR) 
is the dual problem of (PD-STFT-SR). The dual problem of (STFT-SR) is, however, not 
(PD-STFT-SR) as the spaces L^(M^) and Gb(M) are not reflexive. We thus remark that 
the second statement in Theorem 5 concerning strong duality cannot be applied by taking 
(STFT-SR) as the primal problem, that is, with X = Ai{G), Y = L^{G x G), C = Ag, 
(p{u) = ||zz||rpy for all u E Ai{G), and 'ip{c) = X{y}{.c) for all c G T^(G x G). Indeed, the 
function X{y} is not continuous at p, implying that Condition (24) is not met. 

We can now apply Theorem 6 to conclude the following: Take G = M and g{t) = 
;^exp for t E R. Assuming that (PD-STFT-SR) has a solution^, which we de¬ 
note by Co G L°°(M^), the support T = of the measure p to be recovered must satisfy 

|(M*Co)(t£)| = 1, for £ G D, if /i is to be contained in the set of solutions of (STFT-SR). Fur¬ 
thermore, if (PD-STFT-SR) has a solution cq G L°°(M^) such that |x4*co| is not identically 
equal to 1 on M, then every solution of (STFT-SR) is a discrete measure despite the fact 
that (STFT-SR) is an optimization problem over the space of all measures in M(G). This 
can be seen by employing an argument similar to the one in [9, p. 361] as follows: Since the 

condition for the existence of a solution of (PD-STFT-SR) will be given in Corollary 8. 
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window function g is Gaussian and 


Vt e M, {A*gCo){t) = [ co(r, f)g{t - r)e^^*^Mrd/, 

Jr 

both g and A*gCQ can be extended to entire functions which, for simplicity, we also denote by 
g and ^*co, that is. 


E C, g(z) = exp 
'a 


T^Z 

'202 


and 


{A*gCo){z) = 


Co{Tj)g{z-T)e^-^fATdf 


'-fc 


(34) 


It is known that g is holomorphic on C [29, Chap. 10], but it remains to show that so is .4.*Co. 
To this end, let T be an arbitrary closed triangle in C. We have 



(^*co)(^)d2; 


[ ([ [ co{Tj)g{z - T)e^^^f^dTdf] dz 

JdT \J-fc Jr / 

r [ Coir, f )( [ giz - T)e^-^f^dz] drdf 
Jfc Jr \JdT / 

0 , 


(35) 

(36) 


where we used Fubini’s Theorem in (35) and Cauchy’s Theorem [29, Thm. 10.13] in (36). 
To verify the conditions for Fubini’s Theorem, we note that T is bounded, i.e., there exists 
M > 0 such that |; 2 | ^ M for all z E T, which implies that for Al z '■= r + iq E T , r E M., 
and / e [-/c,/c], we have 

\,j(z - !'D 

1 f 7r(r —r)2\ /vrM^h 

where we used |g| ^ M and / E [—fc,fc\ in the last inequality. This yields 
/ /Mr,/)ll9(2-r)e“^*|dTd/ 

= 2/c||co||ioo v^exp^^^^^ exp(2/cM), 
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and hence 


/ [ [ \co{T,f)\\g{z-T)e^^"f^\dTdfdz 

JdT J-fc JR 

< 2L/c llcoll^oo exp(2/eM) < cx), 

where L is the length of dT. By Morera’s Theorem [29, Thm. 10.17] we can therefore conclude 
from (36) that .4.*co is holomorphic on C. We can then dehne the function 


\/z e C, h{z) := 1 - (^*Co)(^)(^;co)( 2 :). 


Since |^gCo| is not identically equal to 1 on M, by assumption, and h is holomorphic on C by 
construction, h is not identically equal to 0. Therefore, thanks to [29, Thm. 10.18] the set 
{z G C: h{z) = 0}, and a fortiori the set {t G M: |(^*Co)(f)| = 1}, are at most countable and 
have no limit points. But since (26) holds, this implies that every solution z/q of (STFT-SR) 
must have discrete support, and therefore, z/q is necessarily a discrete measure. In the case 
G = T, ii g is chosen to be a periodized Gaussian function 


Vt e T, g{t) 





a similar line of reasoning can be applied to establish that every solution of (STFT-SR) is 
discrete provided that (PD-STFT-SR) has at least one solution. 

We have established the existence of a solution of (STFT-SR) and shown that (STFT-SR) 
is the dual problem of (PD-STFT-SR). Now we wish to derive conditions under which the 
measure fi is the unique solution of (STFT-SR). Similarly to [16, Sec. 2.4], we start by 
providing a necessary and sufficient condition for p to be a solution of (STFT-SR), which 
leads to the following theorem. 

Theorem 7 (Optimality conditions). Let fi := ^ A1(G), where Q is a finite or 

eountahly infinite index set and G C\{0}, for all i E Q. Let g G iS(G)\{0}. The measure 
fi to be reeovered is contained in the set of solutions of (STFT-SR) if and only if there exists 
Co G L°°{G X G) such that 

||.A*co||^ ^ 1 and (.A*co)(f^) = |^. (37) 

Proof. Since y = Agfi, it follows by application of [28, Thm. 4.3.6] that p is a solution of 
(STFT-SR) if and only if there exists Cq G L°°{G x G) such that x4*Co G i9||-||rpy (p). Since 
the sub differential of the norm ll-llrpy is given by [16, Prop. 7] 


ViiEM(G), d\\-\\^y{n) = {geGbiG): (h, = II^^IItv > II^IL ^ 1} > 


22 



it follows that /i is a solution of (STFT-SR) if and only if there exists cq G L°^{G x G) 
such that (^*co,/i) = ||/i||Tv ll'^a'^olloo ^ h = condition 

(^*co,/i) = ||/i||Tv equivalent to 

Re i ^ ae{A*gCo){te) [ = • (38) 

I £eo J iefi 

For every £ G ff, we can write 

ae{A*gCo){ti) = \ae\ \ {A*gCo){te)\ , (39) 

where 9i G [0,27r). Thus, (38) is equivalent to 

\ae\ (^1 - cos( 6'£) |(^*co)(t£)| ) = 0. (40) 

e&n 


Since we also have ^ 1, (40) can only be true if |(^*co)(t£)| = 1 and 6^ = 0, both 

for all £ G ff. Therefore, from (39), we can infer that 


Wen, {A*gCo){k) 


ae 

N’ 


which concludes the proof. 


□ 


Thanks to Theorem 7, we can now show the following, which yields a sufficient condition 
for (PD-STFT-SR) to have at least one solution, an assumption made previously to show 
that every solution of (STFT-SR) is discrete. 


Corollary 8. Let fj, := ^ Ai{G), where Q is a finite or countably infinite index 

set and G C\{0}, for all £ e n. Let g G 5(G)\{0}. ///i is contained in the set of solutions 
o/(STFT-SR), then (PD-STFT-SR) has at least one solution. 

Proof. If the measure g, to be recovered is a solution of (STFT-SR), then by Theorem 7, 
there must exist Cq G L°°{G x G) such that (37) holds. Since ^ 1, Cq is feasible for 

(PD-STFT-SR). Moreover, this cq satishes 

{co,y) = {co,Agfj,) = {A*gCo,fi) = Re |^ {A*gCo){t)fj,{t) 

= Re<^ ^afiA*gCo){tfii I =^\ae\ = ||/i||Tv> 

I e&Q J £€Q 

which shows that the supremum in (25) is attained for cq. Therefore, cq is a solution of 
(PD-STFT-SR). □ 
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Theorem 7 provides conditions on yU to be a solution of (STFT-SR). However, we need 
more, namely, conditions on /i to be the unique solution of (STFT-SR). Such conditions are 
given in the following theorem, which is a straightforward adaptation of [3, App. A], 

Theorem 9 (Uniqueness). Let fx ;= ^ M.{G), where VL is a finite or countably 

infinite index set and ag G C\{0}, for all £ G fl. Let g G iS(G)\{0}. If for every sequence 
e = of unit magnitude complex numbers, there exists a function Cq G L°°{G x G) 

obeying 


'il G ff, (.4.*co)(U) — Si (41) 

VfGG\T, |(^*co)(t)| <1, (42) 

where T := then p, is the unique solution of (STFT-SR). 

Proof. Assume that there exists a Cq G L°°{G x G) such that (41) and (42) are satisfied for 
all sequences of unit magnitude complex numbers. Let z/q G M.{G) be a solution 

of (STFT-SR) and set h := p — uq. Suppose, by way of contradiction, that h 7 ^ 0. The 
Lebesgue decomposition [29, Thm. 6.10] of h relative to the positive and a-hnite measure 
= 'Eeen K is given by 

h = hj' -|- hj'c, 

where hr is a discrete measure of the form hr '■= hiSt^ (see Part (b) of [29, Thm. 6.10]) 
and hT<= is a measure supported on := G\T (see Part (a) of [29, Thm. 6.10]). Using the 
polar decomposition [29, Thm. 6.12] of hx, we hnd that there exists a measurable function g 
such that \g(t) \ = 1, for all t eT, and 

[ ^{t)hT{t) = [ (p{t)g{t)\hT\{t) (43) 

Jg Jg 

for all (p G Gb{G). By assumption, there exists a cq G L°°{G x G) such that 

G H, (.A*co)(U) = ?7(U), (44) 

VtGG\T, |(^;co)(t)|<l. (45) 
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We thus have 


(^*Co, — (^*Co, h) — (^*Co, hx) 

= {cQ,Agh) -Re| /” (^*co)(t)hT(t) 


IG 


= {co^Agii - AgPo) - Re |y (^;co)(t)77(t)|hr|(t) 

= - Re i ^ \hf,\ {A*gCQ){t^)ri{t^) I = - ^\h(,\ \r]{te)\ 

I ten J 

= ~ \hi\ = — ||hr||Tv ’ 


£en 


ten 


(46) 

(47) 

(48) 


where (46) makes use of (43), (47) follows from (44) and y = AgH = AgUo, which holds 
because z/q was assumed to be a solution of (STFT-SR), and (48) is by \ri(t)\ = 1, for all 
t & T . It therefore follows that 


Ik^rl TV “ \{'^*gCo, hTc)\ — 

Re| / {A*gCo){t)hTc{t)\ 


[ {A*gCo){t)hTc{t) 


Ug ) 


Jg 

< jji^*gCo){t)\ \h 

te (^) 




< \hTA (^) — II II TV ’ 


(49) 

(50) 


where (49) is by the triangle inequality and the strict inequality in (50) is due to (45) and 
h 7 ^ 0 and hence ^ 0 (as = 0 would imply hr = 0 and thus h = 0, because of (49)). 
We then obtain 


II^oIItv ~ 11^ ~ ^IItv ~ 11^ ~ ^tIItv II^t=|Itv 

^ Ilh-llTV ~ II^'tIItV II^'T^IItV 
> IIfIItv ’ 

where (51) is a consequence of y being supported on T and hence the supports of y — hx and 
h-rc being disjoint, (52) follows from the reverse triangle inequality and (53) is a consequence 
of (50). This contradicts the assumption that uq is a solution of (STFT-SR), which allows 
us to conclude that h = 0 and hence z/q = y, thus completing the proof. □ 

We have now developed a fnll theory of snper-resolntion from STFT measurements for 
window functions g from the Schwartz-Bruhat space that are extendable to entire functions. 
It remains, however, to connect the nniqneness conditions (41) and (42) in Theorem 9, which 
amonnt to constrained interpolation problems, to the minimnm spacing condition A > l//c, 
as a nn ounced earlier in the paper. This will be accomplished throngh a specihc choice for g, 
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namely we take it to be a Gaussian for G = M and a periodized Gaussian for G = T. Note 
that the choice of the window width is only for reasons of specihcity of the results. Similar 
results can be derived for other Schwartz-Bruhat window functions or for Gaussian window 
functions of smaller widths: we only need to adapt the computational part of our proof, but 
the way of reasoning remains the same. 


Theorem 10 (Gonditions for exact recovery in the case G = M). Let G = M and 


Vt G M, 






where (r ^■ Let fi := ^ A1(M), where LL is a finite or countably infinite index 

set and G C\{0}, for all i E fl. If the minimum distance 


A := min \ti - tm\ 

satisfies A > 1/fc, then the conditions of Theorem 9 are met, and hence fi is the unique 
solution o/(STFT-SR). 

The proof architecture of Theorem 10 is inspired by [3, Sec. 2, pp. 15-27]. There are, 
however, important differences arising from the fact that we deal with STFT measurements 
as opposed to pure Fourier measurements. Specihcally, in the case of pure Fourier measure¬ 
ments, the interpolation function v4.*Co must be a Paley-Wiener function [29, Thm. 19.3], 
while here .4.*co is clearly not band-limited and can therefore have better time-localization. 
We believe that this allows the minimum spacing A to be smaller than in the case of pure 
Fourier measurements. 

For G = T, we have /c = IF G N. Moreover, the set G has to be hnite, as T is compact 
and G indexes T which is closed and discrete. Recovering the measure fi therefore reduces to 
the recovery of the hnite set of points C [0,1), L := |G|, and the associated weights 

{orjteo- 

Theorem 11 (Gonditions for exact recovery in the case G = T). Let G = T and 

Vt G T, ^(t) ;= ^ exp (-j , 

where a := 4 (x+i/ 2 ) • h' G A1(T) with ae G C\{0}, for all i E {1, 2,..., L}. 

If the minimum wrap-around distance 


A := min min \te — tm + n 
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satisfies A > 1/{K + 1/2), then the conditions of Theorem 9 are met, and hence /i is the 
unique solution o/(STFT-SR). 

6 A recovery algorithm for G = T 

We next provide an explicit recovery algorithm for the case G = T. Specihcally, we show that 
if fi is the unique solution of (STFT-SR) (which is the case provided that A > 1/{K + 1/2) 
and a = l/(4(iF+l/2))), then an approximation of the correct solution fi can be recovered by 
solving the convex programming problem (STFT-SRat), N eN, dehned below. The predual 
of (STFT-SRtv), unlike that of (STFT-SR), is equivalent to a hnite-dimensional problem, 
which can be solved numerically. The justification for this procedure is given in Proposition 
3, which shows that the sequence of solutions un of (STFT-SRat) converges in the weak-* 
sense to /i as TV —)■ cx). 

We hrst note that the 1-periodic window function g can be expanded into a Fourier series 
according to 

Vt e M, g{t) = ^ exp ^ 

nGZ ^ ^ / nGZ 

with g^ = v^^exp (—27r(T^n^), n G Z. The corresponding STFT measurements of fi are 
given by 

WeT,Wke{-K,...,K}, y{T,k) = {Vgfi){T,k)= 

Jo 

L 

l=\ n<^'L 

where yk,n '■= J2f=i is the nth Fourier series coefficient of r h-)■ y(T, k). Using 

Parseval’s theorem, the objective function in (PD-STFT-SR) can be rewritten as 

(c, y) = Re ^k,nykA ^ , 

Ifc=-Xnez j 

where Ck,n denotes the nth Fourier series coefficient of the function r h-)■ c{T,k), for k G 
{—A,..., iV}, and c G L°°(T x Z) is the optimization variable of (PD-STFT-SR). For 
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c G L°°(T X Z), the function A*c is given by 


K „i 


VteT, {A*gC){t)= ^ [ c{T,k)9{t 

k=-K-J^ 

= Y. rc(r,A;) 

k=-K'^^ VneZ / 

EE gn( f c(r, A;)e-2™Mr') 

h — —TrnCW \J0 / 


k=—K nSZ 
K 

EE 

k=—K nGZ 


9n^k,n^ 


27 vi{k-\-n)t 


(55) 


Since infinitely many coefficients Cfc^„ are involved in the expression (55) of the feasible 
set {c G L°°(T X Z): ^ 1} of (PD-STFT-SR) is infinite-dimensional. We now 

approximate this infinite-dimensional problem by retaining only 2N + 1 coefficients in the 
Fourier series expansion of g, that is, we replace g by 

N 

Vt G T, g^{t) = Y 

n=—N 


where N is chosen large enough for the coefficients gn = v^^exp (—27r(T^n^), |n| ^ iV, to be 
“small”. The problem (STFT-SRat) is now defined as 


(STFT-SRtv) minirnize ||zz||rpy subject to y = Ag^u. 

As Theorems 4 and 6 hold for every g G iS(T) and we have g^ G 5(T), they can be ap¬ 
plied to conclude that (STFT-SRat) always has a solution. The Fenchel predual problem of 
(STFT-SRat) is given by 

(PD-STFT-SRat) maximize {c,y) subject to \\AX ell ^1, 

where the measurements now become y{T,k) = (‘^siv„,^w™)('r, ^)t for r G T and /c G Z. 
Moreover, strong duality holds, implying that 

min I llz/ll^y : Ag^^u = |/, z/ G Af(T)| = sup | {c,y) : ^ l,c G L“(T x Z)|. (56) 

Now, we make the problem (PD-STFT-SRat) explicit. The objective of (PD-STFT-SRtv) can 



be written as 


:,y) = {C,Y) = Re\ E 


^k^nVk ,n 


k=-Kn=-N 


where C ;= {ck,n)\k\<^K,\n\^N and Y := (|/fc,n)|fc|=£ii',|nKAf- The function A~^c is a trigonometric 
polynomial (hereafter referred to as dual polynomial), entirely characterized by the hnite set 
of coefficients {cfc,n}|fe|^ic,|n|^Af, and expressed as 


VieT. (^yc)(()= 5: 5: 


K N K+N 

E E ■ 

k=—K n=—N m=—{K+N) 


where 


E 9mCm-n,n, with n^in ^ = max{- Y, 771 - 


J^max := min{Y, m + K}. 

The problem (PD-STFT-SRat) thus takes on the form 


maximize (C,Y) subiect to \\m=-(K+N) 

C(zC(2K+1)x(2N+1) II I loo 

9nCm—n,n- 

^ "^min 

Now, as in [3, Sec. 4, pp. 31-36], we make use of the following theorem to recast 
(PD-STFT-SRat) as a semi-dehnite program. 

Theorem 12 ([30, Cor. 4.25]). Let M eN and let P be a trigonometric polynomial 


Vt e T, P{t) = P- 


Then, ||T’||^ ^ 1 */ and only if there exists a Hermitian matrix Q G (^(2M+l)x(2M+l) 
that 


Q p 

1 


0 and Qm, 


m+i — 


1 , £ = 0 

0, £g{1,...,2M}, 


where p G is the column vector whose kth element is pk-M-i for /c G {1,..., 2M + 1}. 
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From (57) and Theorem 12 it then follows that (PD-STFT-SRat) is equivalent to 


maximize (Y, C) subject to 

CeC(2^f+l)x{2]V+l) 

QgC*^'xM' 

where M' := 2{K + N) + 1, x is the column vector whose kih element is xu-k-n-i for 
k E {1,..., M'}. We can now solve (PD-STFT-SRtv) and recover the corresponding measure 
by polynomial root finding as done in [3, Sec. 4, pp. 31-36]. 

The following result shows that the sequence of solutions of (STFT-SRtv) converges in 
the weak-* sense to /r as iV —)■ cx), provided that fi is the unique solution of (STFT-SR). 

Proposition 3. Let fi := ^ W1(T) with ae E C \ {0}, for all ^ E {1, 2,..., L}. 

If is the unique solution o/(STFT-SR), then every sequence {zzjvItvsn such that for every 
N eH, vn is in the set of solutions of (STFT-SRat) converges in the weak-* sense to jj,. 


Qn^m—r 

lin 

P 0 


Q X 

x^ 1 


( 58 ) 


M'-e 

Qk,k+i = ^0/, 

k=0 


Proof. Let N eN. By (56) we have 

ll^ivIlTv = sup |(c,|/) : ^ ^ ^ ^)} • 

For c G L°°(T x Z), thanks to (57), {A~^c)(t) depends only on the Fourier series coefficients 
Ck,m n G {—N ,..., N}, of the functions r h-)■ c(r, k), for k E {—R",..., 77}. The problem 
(PD-STFT-SRat) is therefore equivalent to 


maximize (c, y) subject to ||x4,*cl| ^ 1, 
cepTv(TxZ) ' •’ II S Iloo ^ ’ 


where 


'PAr(T X Z) := < (r, k) h-)- 


N 


0, otherwise 


-N€n€N 


But since 

(c.v) ■ ||xljc||^ < l,ce X Z)| C |{c, 9 > : ||xl’c||^ < 1, c e'Pw+i(T 

c{(c,j/): llxiyll^^ l,ceL“(TxZ)}, 


X 
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we have that 


as well as 


\^n\ 


TV 


= mm 



II ^ 1, c 

11 oo ^ 

G 7^A''('ir X ^) 



■ IMy 

II ^ 1, C 

11 OO ^ 

£ 'Pn+iO^ X Z)| 

— Ik'V+ll TV 

(59) 

{(c,!/) : 


^ 1, c G P]\fiT X 

Z)} 


{{c,y} ■ 

IM;c|L 

^ l,c G L“(T X 

Z)} 


{| ^1 TV 

: Agiy = y,iy e A1(T)} = || 

FI TV • 

(60) 


From (60) it follows that the sequence {z/ArjArg^ is bounded. Therefore, by application of [23, 
Cor. 3.30], there exists a subsequence {z^Ar^^jmeN that converges in weak-* topology to a 
measure Uoo G A^(T). From (59) and (60) if follows that {Hz^vIlTvi-^^eN is convergent. Thanks 
to [23, Thm. 3.13 (hi)] and (60), we then get 


1 ^ 00 11XV ^ lim llz^Vmllxv ^ 

m^oo 


TV 


Next, we show that 

lim (AgUNjir, k) = |/(r, k). 

m^oo 

Let r G T and k G {—K,... ,K}. For m G N, we have 

\y{T,k) - {AgUNj{r,k)\ = \{Ag^^uNj{T,k) - {AgUNj{T,k)\ 

r 

—27Tikt 


9Njt - - git- r)e-^-^V^„(t) 


{dN^it - r) - g{t - r))e 


< / lf/v™(t-r)-^(t-r)l |z/Ar,„| (t) 

Jt 

Iljiv„ - jIL IIi'n.JItv 


(61) 


2 2 \ 
' n ] 


From 

CXD 

WdN^ - ^IL ^ 2 ^ v^exp(-27ra2 

n=Afm+l 

it follows that lim^^oo \\gNm “S'lL = 0- Since {Hz/w™ ||tv}"i 6N is also convergent, (61) con¬ 
verges to 0 as m —)■ cx). It therefore follows that 


lim [AgVNjiT, k) = |/(r, k). 
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But since {z^Ar„}meN converges to z/qo in the weak-* topology and 

Jj 

we have, by dehnition of weak-* convergence, 

lim {AgUNm){'T,k) = {AgUoo){T,k)- 

m^oo 

We therefore get y{T, k) = {AgUcAji^, k), for all r G T and k G {—K, ..., K}, which shows 
that z/qo is feasible for the problem (STFT-SR). As y is the unique solution of (STFT-SR), 
we have ||/i||rpy ^ ||^oo||tv’ which combined with ||z/oo||tv ^ IIfIItv l^nds to ||z/oo||tv ~ II/^IItv’ 
and hence shows that z/qo is a solution of (STFT-SR). However, by assumption /i is the unique 
solution of (STFT-SR). We can therefore conclude that not only ||z/oo||tv ~ II/^IItv 
T^oo = y- Since there was nothing specihc about the accumulation point z/qo, we can apply the 
same line of arguments to every accumulation point of the sequence {z/Ar}AreN, and therefore 
conclude that every accumulation point of {z/Ar}ArgN must equal /i. In summary, we have 
shown that /i is the unique accumulation point—in the weak-* topology—of the sequence 
But since ||z/Ar||rpy ^ H/^IItv (60), the sequence {z/atItvsn is contained in the 
closed centered ball of radius ||/i||TV) which, by the Banach-Alaoglu Theorem, is compact in 
the weak-* topology. We next show that this implies weak-* convergence of {z/Ar}ArgN to y. 
Suppose by way of contradiction that {z/ArjArg^ does not converge—in the weak-* topology— 
to /i. Then, there exists an £ > 0 such that for inhnitely many N we have \{un — 
for all if G C(T). Therefore, we can hnd a subsequence of {z/Ar}AfeN which, by compactness— 
in the weak-* topology—of the closed centered ball of radius ||yn||Tv; accumulation 

point different from /i in the weak-* topology. This constitutes a contradiction and thereby 
hnishes the proof. □ 

Fix A^ G N. If (PD-STFT-SRtv) has at least one solution c G L°°(T x Z) such that 
|^~^c| is not identically equal to 1, it follows by (26) that, as a consequence of being 
a trigonometric polynomial, every solution z/at of (STFT-SRat) is discrete. Unfortunately, 
the weak-* convergence alone of the sequence {z/Af}AreN of discrete measures to the discrete 
measure /i, as guaranteed by Proposition 3, does not imply, in general, that each element of 
supp(z/Ar) converges to an element of supp(/i). What we can show, however, is that for small 
enough e > 0, one can hnd an > 0 such that for all N ^ iVA, each set [ti — e,t£ + e], 
i E {1, 2,..., L}, contains at least one point of supp(z/Ar), and in addition, most of the “energy 
(in TV norm)” of pn is contained in := IJnez U^=i [te + n — e,te + n + e]. This is formalized 
as follows. 

Proposition 4. Let /i = £ A4(T) with ae E C\ {0}, for all i E {1, 2,..., L}. 
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Assume that the wrap-around distanee 

A := min min \t£ — tm + n\ >0. 

If fi is the unique solution 0 /(STFT-SR), then every sequenee {z^atItvsn such that for every 
iV G N, z/tv is contained in the set of solutions of (STFT-SRat), satisfies 

We E (0, A/4], 3Ms > 0, \/N ^ Wi E {1,2,... ,L}, supp(z/Ar)n[tf —e,t£ + e] 7 ^ 0. (62) 


Moreover, setting 


We E (0, A/4], Tg [tf — e + n,+ e + n] 

nGZ £=1 

T/:=T\T„ 


and defining i'N,Te)^N,T^ G A1(T) according to 


WB E B{T), vn,tAB) ■■= MTe n B) 

:= n B), 


we have 


We E (0, A/4], 3M, >0, WN ^ 


I 11 rpY ^ £ IIPIITV 

|^Af,rJlTV 


(63) 


TV 


Some remarks are in order before we prove Proposition 4. An obvious consequence of 
(62) is that the number of atoms of z/jv is larger than (or equal to) the number L of atoms of 
fx. Moreover, for all N ^ and t^^'^ E supp(z/Ar) fl T/, we have 




^ Piv,r-||Tv 

^ £ IIfIItv ’ 


(64) 

(65) 


where (64) follows from e Tf and (65) is by ||z/jv,r=||Tv = Ezesupp(^^)nr= kMR({f})|- 
This means that for sufficiently large N, the weights of un attached to the spikes located 
outside Tg are smaller than e ||yu||rpy, provided that e G (0,A/4]. Note, however, that each 
of the sets [t£ — e,ti -\- e], i E (1,2,..., L}, may contain more than one atom of z/^r. Fig. 2 
illustrates the statement in Proposition 4. 

Proof. We start by establishing (62). To this end, let e G (0,A/4], i E (1,2, ...,L}, and 
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t3-£ ^3 *3+5 


tl-E tl <l+£ 


* 2—6 t 2 *2 +£ 


Figure 2: Illustration of Proposition 4 with L = 3. The union of intervals depicted in red 
represents Tg on the fundamental interval [0,1). The spikes in red represent fi. Proposition 
4 guarantees that for N ^ Mg, each interval + e\ contains at least one atom of 

The spikes in blue correspond to for N ^ Mg. The weights of the atoms of z/at outside Tg 
is guaranteed to be smaller than e ||/i||rpy. 
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define G C(T) as 


Vt e T, 

where : M —)■ M is given by 




Wei 


nGZ 


t — te — n 


Vt e M, 



\t\^l 

otherwise. 


( 66 ) 


As the minimum wrap-around distance between the te is A, by assumption, and e ^ A/4, 
we have 

L 

amlpeAtm) = at'lpeA'^d = W\ ■ 

m=l 

As {z^vjAfeN is, by assumption, a sequence of solutions of (STFT-SRat) and fx is the unique 
solution of (STFT-SR), it follows from Proposition 3 that {z^atItvsn converges in the weak-* 
sense to p., i.e., for every ip G C(T), it holds that 


V?7 > 0, 3M^^ri >0, VA ^ \{<P,i^n) - (9^,h)| < V- (67) 

Since tp is continuous, by construction, G C(T), and setting rj = lap /2 > 0 and ip = ip^^i 
in (67) implies that there exists an Mi;^£ > 0 such that for all N ^ we have 

\{'4’e,i,^N) - (^M,h)l = - la^ll ^ \a£\/2. (68) 


Now, set 

Me := max Me/, 

let N ^ Me, and assume, by way of contradiction, that there exists an i' G {1,2,...,L} 
such that supp(z/Ar) fl [t// — e,t£/ + e] =0. Then, as p’s/' = 0 on [0, l)\[t/' — e,t£/ -1- e], we 
have {ipe/'TixN) = 0, which by (68) would imply |a//| ^ |a/'|/2 and thereby contradict our 
assumption a// ^ 0. We can therefore conclude that for N ^ Me, for all 7 G {1, 2,..., L}, we 
have supp(z/Ar) fl [f/ — e,t£ + e] ^ 0. This completes the proof of (62). 

We proceed to establishing (63). For e G (0, A/4], dehne tp^ G C(T) as 


L 

VtGT, ^/>e(t) := ^^£,£(t), 

£=1 


where tpe^i'. M —)■ M, 7 G {1,2,...,L}, is as in (66). As for all 7 G {1,2,...,L}, ip^^£ is 
supported on — e + n,t£ + e + n], tp^ is supported on T^. Moreover, since e ^ A/4, 
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by asssumption, the ipe/ have disjoint supports, and hence 


L 

Vt e T, \pje{t)\ = ^\i^sAt)\ 

e=i 


(69) 


We also have 


L L 

(Aei h) ~ ^ ^ h) = ^ ^ = Ilh-llTV ‘ 

e=l 1=1 

Applying (67) with (p = and rj = £||/i||rpy > 0, we can conclude that there exists an > 0 
such that for all N ^ we have 


^ II/^IItV ^ \{'^e,^N) - (^£,h)l 

= li'P’e, ^N,T^) — IlhliTvl 

^ 11/^11 TV ~ {'^ei ^N,Ts) ) 

where (70) is a consequence of tp^ being supported on T^. It therefore follows that for all 
^ Mg, we have 

(1 - e) ll/illxv ^ ^ / \'P’e{t)\\l'N,TMt) ^ [ WN,TMt) = \An,T,\\tY ^ 

Jo Jo 

where the last inequality follows from (69). Furthermore, we have for all N ^ Mg, 

^ II^AfllxV ~ (1 ~ Ilh-llTV 

^ Ilh-llTv ~ ~ IIfIItv ~ ^ IIfIItv ’ A A 

where (72) is a consequence of the supports of and i'n,t^ being disjoint, (73) follows 
from (71), and (74) is by Ur'ArllTv ^ II/^IItv’ iV G N, as established in (60). □ 


7 Numerical results 

For our numerical results we consider the case G = T, i.e., recovery of the discrete complex 
measure /i = ^ A1(T) from the measurements |/(r, k) = (^g^z/Ar)(r, k). We solve 

the predual problem (PD-STFT-SRat) with N = 25 and iV = 50 by applying the convex 
solver cvx to the formulation of (PD-STFT-SRat) as given in (58). 

To assess recovery performance, we run 1500 trials as follows. For each A, we construct 
a discrete complex measure /i supported on the set T = {G}f=o with S = [1/(2A)J and 
tg = 2iA + rg, where rg is chosen uniformly at random in [0, A]. The minimum wrap-around 
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Figure 3: Success rate for support recovery from pure Fourier measurements and from STFT 
measurements with fc = 50 and = 50. 

distance between the points in T is therefore guaranteed to be greater than or equal to A. The 
complex weights a£ are obtained by choosing their real and imaginary parts independently 
and uniformly at random in [0,1000]. We declare success if the reconstructed measure fi 
has support T = satisfying ||T — T||£ 2 /||T ||£2 ^ 10“^. The corresponding results 

are depicted in Fig. 3. We observe that both recovery from STFT measurements and from 
pure Fourier measurements actually work beyond their respective thresholds A > 1/fc and 
A > 2/fc, thus suggesting that neither of the thresholds is sharp. We also observe a factor- 
of-two improvement in the case of recovery from STFT measurements relative to recovery 
from pure Fourier measurements, suggesting that the improvement in the recovery threshold 
A > l//c for STFT measurements relative to A > 2//c for pure Fourier measurements is 
due to the recovery problem itself. Specifically, in the STFT case, we perform windowing 
and the STFT measurements are highly redundant. To see this note that in the pure Fourier 
case, the measurements consist of the vector y = {yk}k=-K ^ (^2K+l 

as defined in (6), 

whereas the STFT measurements here are characterized by the {2N + 1) x {2K + 1) entries 
of the matrix Y G containing the 2N + 1 Fourier series coefficients of the 

2K + 1 functions r h->■ y{T, k), k E {—K,... ,K}. The increased number of measurements 
in the STFT case leads to larger optimization problem sizes and hence entails increased 
computational complexity relative to the pure Fourier case. 

In Fig. 4, we compare the dual polynomial for pure Fourier and for STFT measurements 
in a situation where recovery in the former case fails and where it succeeds in the latter.^ We 

^In the case of pure Fourier measurements, we follow the recovery procedure as described in [3] and solve 
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can see that in both cases, the magnitude of the dual polynomial equals 1 for all te, £ G hi. 
However, in the case of pure Fourier measurements, the magnitude of the dual polynomial 
also takes on the value 1 at eight additional locations not belonging to the support set 
T. For example, we can see in the top plot in Fig. 4 that spikes are detected for t = 
0.1791 and t = 0.2172, while these two points do not belong to the support set of the 
original measure. In contrast, in the case of STFT measurements, the locations where the 
magnitude of the polynomial takes on the value 1 approximate the points of the support set 
of the original measure (represented by circles in Fig. 4) well, namely with a relative error of 
||f-T||,V||T||,. = 2-10-7. 

A Proof of Theorem 10 

Let e = be a sequence of complex unit-magnitude numbers. The goal is to construct 

a function cq G such that (^*co)(t£) = for all £ G fl, and |(^*co)(t)| < 1 for all 

t G M \ T, where T = is the support set of p. Inspired by [3], we take cq to be of the 

form 

V(r, /) G co(r, ^ {ai9{k - + (it g{ti - , (75) 

where (3e G C, for all i E fl, and we proceed as follows: 

1. We first verify that a := {aejeen and {3 := both in £°°(r2) implies Cq G L°°(M^). 

2. Next, we show that one can find a,/3 E such that the interpolation conditions 

{A*Co){te) = ££, for all i E il, are satished and |-d,*Co| has a local extremum at every te, 

£ eQ. 

3. Then, we verify, with a,/3 E £°°{fl) chosen as in 2., that the magnitude of A*gCo is 

indeed strictly smaller than 1 outside the support set T = of p. This will be 

accomplished in two stages. First, we show that |^gCo| is strictly smaller than 1 “away” 
from each point ti, £ E Q, specihcally, on 7/c , t£ + ^]. We then complete 

the proof by establishing that |-d,*Co| is strictly concave on each set [t^ — ^,te + 

£ E Vt, which, combined with the fact that |(^*co)(f£)| = 1, for every £ E Vt, implies 
that |^gCo| is also strictly smaller than 1 on each of these sets. 

The main conceptual components in our proof are due to Candes and Fernandez-Granda 
[3]. Although [3] considers recovery of measures on T only and from pure Fourier measure¬ 
ments, we can still borrow technical ingredients from the proof of [3, Thm. 1.2]. However, 

(D-SR) defined in Section 3 using the convex solver cvx. 
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Figure 4: Example of reconstruction of a discrete measure fi with |supp(p)| = 54 and A = 
0.011880 from both pure Fourier and STFT measurements with = 50 and N = 25. The 
curves represent the dual polynomials and the circles represent the elements in the original 
measure’s support set. 
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the different nature of the measurements and, in particular, the case G = M, pose additional 
technical challenges relative to the proof of [3, Thm. 1.2], Specihcally, the sum corresponding 
to (75) in [3] is always hnite, whereas here it can be infinite, which presents us with deli¬ 
cate convergence issues that need to be addressed properly. Further fundamental differences 
between the proof in [3] for the pure Fourier case and our proof stem from the choice of 
the interpolation kernel, which here is given by t h-)■ i?(t) sinc(27r/ct). Specihcally, we do 
not have to impose a bandwidth constraint on the interpolation kernel. For pure Fourier 
measurements, on the other hand, the interpolation kernel has to be band-limited to [—/c, fc] 
(Candes and Fernandez-Granda [3] use the square of the Fejer kernel which offers a good 
trade-off between localization in time and frequency). As already mentioned in the main 
body, this leads to a factor-of-two improvement in the minimum spacing condition for STFT 
measurements over pure Fourier measurements. Note, however, that STFT measurements, 
owing to their redundancy, provide more information than pure Fourier measurements. We 
hnally note that our proof also borrows a number of technical results from [17]. 

A.l a,(3 ^ implies cq G L°°(R^) 

Let a := G ^“(G) and {3 := G t'°°(r2). Take (r,/) G and dehne IqRi G G 

to be the indices of the points in T, that are closest and second closest, respectively, to r, 
that is, 

£0 := arg min IG — t| and := arg min |G — r| . 

For brevity of exposition, we detail the case ^ r ^ only, the cases ^ r ^ ^ 

Go ^ r, and r ^ Go ^ dealt with similarly. For all £ G G7 := {m E fl: tm < Go}) 

it then holds that 

T - ti ^ ti^ - ti ^ A, (76) 

and for all i G := {m Gi), we have 

te - T ^ ti - ^ A. (77) 
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Hence, we get the following: 


|co(r,/)| ^ ^ I \\a\\i^^\gite - t)\ + \\/3\\^^'^\g'{te-T)\] 

\ ien een / 

= ^ |j(*4 - 01 + |J(0 - 01 + \g'(t,, - T)\ 

+ !!f|^l9'(ff.-OI + ^|f^ E l9(«f-OI + Jlf|^ E l9'(«'-OI 

£eo\{4,£i} £eO\{4,£i} 

< Ihlt iijii^ + MJt ||g'||_^ + !!^ E Isfe. - *')! 

teOr 

+ Ihjf E Is!*' ^ E Is'!*'. - <')l + ^ E !»'(«' - 0)1. 


(78) 


(79) 

(80) 


where the step from (78) to (79)-(80) follows from g,g' G Cfe(M), (76), (77), and the fact that 
and \g'\ are both symmetric and non-increasing on [A, cx)), which is by the assumption 
A > 4a. Note that we eliminated the dependence of the upper bound in (79)-(80) on (r, /). 
It remains to establish that every sum in the upper bound (79)-(80) is hnite. The minimum 
separation between pairs of points of T = is A, by assumption. Consequently, since 

and \g'\ are both symmetric and non-increasing on [A, cx)), the sums in (79) and (80) take 
on their maxima when the points £ G ff, are equi-spaced on M with spacing A, i.e., when 


{tio — ^ G ff.,- } C |nA; n G N \ {0}} 
{tt — tq : f' G } C |nA: n G N \ {0}}. 


It therefore follows that 


a 


|co(r,/)|^^ 

Jc 


where 


+ hfc lls'lL + 


fc 


fc 


Elj(nA)| + f^El9'(nA)|, (81) 


n=l 


n=l 


n=l 

oo 

Ei9'("A)i 


Eexp 


71=1 


Tin 


2^2 


2a2 


TT 


n=l 


cr^x/a 


n A exp 

n=l 


< OO 


Trn^A^ 
2a2 


< cx. 


which establishes that Cq G L 


oom2'i 
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A.2 Existence of G such that {A*gCo){t£) = S£ and |A*co| 

has a local extremum at every G U 

Using (75) in (18), we get 

VteM, {A*gCo){t)=[ /'co(r,/)^(t-r)e^”^Mrd/ 

J-fc Jm. 

^ f I £rYl - T)e-^’“f“ )g(t - T)e="-"dTd/ (82) 

+ft g'(t, -T)g(t- r)dTj 1'° (84) 

= ^ («£ R{t - te) sinc(27r/c(t - te)) +13^ B!{ti - t) sinc(27r/c(t - ti)) ^, (85) 


where we set 


■.=u(t-ti) 


Vt G M, u{t) := R{t) sinc(27r/ct) 




Vt G M, v(t) := -R'(—t) sinc(27r/ct) = ^-^i?(t) sine(27r/ct) 


Vt G M, i?(t) = exp 

as defined in (15). The conditions for Fubini’s Theorem, applied in the step from (82)-(83) 
to (84), can be verified as follows: 


/ / \c^,9{k - T)g{t - r)e2-/h-*d| drd/ 

J-fo -'K 

[ \9{te-T)g{t-T)\dT] ( 7 ^ [ 
fan V./K / \^Jc J-f, 


Y, l«£l R{t - U) ^ 2 ||a||,^ ||i?i|^ + i?(nA) 
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and 


/ / \Pe9'{te-r)g{t-T)e‘^^"^^^ *^^|drd/ 

e&n ''-f^ 


£en 


/ \9{te-T)g{t-T)\dT ( df 


1 


2fcJ- 


fc 


^ Y, m m -t)^2 ni^iu + R{n^) ), 

£gr 2 V n=l 


where we used 


VteM, y'^(r) |^'(t + r)|dr = ^exp|^-|^ 


a-w'-'is- 

< + |fi'(()| =: R(t) 


(89) 


and the fact that R is bounded, symmetric, and non-decreasing on [A, cx)) as a consequence 
of A > da. The upper bounds in (88) and (89) are both hnite as the series Yhn^i RiP'^) and 
R{nA) converge. 

We have shown in Lemma 3 that for cq G L“(M^), the function A*gCo is in C'fe(M). With 
Co taken as in (75), ^*Co is not only in Cb(M), but also differentiable, as we show next. We 
start by noting that the functions u and v dehned in (86) and (87) are differentiable on M, 
and their derivatives are given by 


Vf G M, u' it) = R!{t) sinc(27r/ct) -|- 27r fcR{t) sinc'(27r/ct) 

Vf G M, v'{t) = —R"{—t) sinc(27r/ct) -|- 27rfcR'{—t) sinc'(27r/ct) 

smc(27r/ct) + —^smc (27r/ct) 


2a‘^ da^ 


a^ 


R(t). 


Then, using 


Vt G M \ {0}, |sinc(t)| ^ — and |sinc'(t)| = 

rl 


cos(t) sin(t) 




1 1 

< ^ (90) 

i«i iir 


we obtain the following upper bounds on u' and v'\ 
VtGM\{0}, \u'{t)\^ 

VtGM\{0}, |n'(t)| ^ 


1 1 
+ ^ + 


da^fc ' \t\ ' 27rfc\t\^ 
1 


R{t) =: U{t) 


+ A + TTI R(t) =-. v(t). 


(91) 

(92) 


20-2 /c | t | 2 o -2 80-4 /c 

Next, we establish that Yhem ~ ^e) + — ti)) converges uniformly on every com- 


43 



pact set [—r,r], r > 0, so that we can apply [31, Thm. V.2.14] to show that the series in (85) 
can be differentiated term by term. For r > 0, we have 


sup \aiu'(t - ti) + I3ev\t - ^ Ikll^oo sup \u'{t - ti) \ + \\/3\\^^ sup \v'(t - ti)\ 

g(,^te[-r,r] \ te[-r,r] te[-r,r] 


lalhoc 1 Z-, sup \u'{t — t£)\ + i:. sup \u{t — t£)\ + sup \u'{t — t£)\ 




+ ll/^ll<;- I X] sup \v'{t-t£)\+ sup \v'{t-te)\+ sup \v'{t-t£)\\, 

. eeUr *£[-»•,d ie[-r,r] te[-r,r] 




ieflr 


where we dehned the sets := {i eQ: t£ E [—r,r]}, := {i E fl: te > r}, and hi” := 

{i E Q: te < —r}. The functions U and V are both positive and symmetric, U is non¬ 
increasing on (0, cx)), and V is non-increasing on (0, ex) as 


Vte(0,x), V\t)=(^- 


+ 


TT 


2a^U^ 8aV, 


R{t) + 


1 TT 

+ 


1 TT 

+ + 




2a‘^fct 2a‘^ 
t + 




2a‘^fcR Sa'^fc \2a‘^ / IGcr^fc 


^ 8(j4/c/ 

R{t) < 0. 


It therefore follows that 


W 


sup \u\t-t£)\ ^ ^U{r-tf), t£>r 

t£[—r,r] 


and 


{ ll^'1loo> t£E[-r,r] 

V{r-t£), t£>r 

V{-r-te) t£<-r. 


Since the support set T = is closed and uniformly discrete, by assumption, and [—r, r] 

is compact, the set T fl [—r, r], and thereby the index set contains a hnite number of 
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elements, say Lr G N. We thus have the following 


/ ^ sup \atu{t - te) + (3ev\t - te)\ ^ ||q;||^oo -hr IIm'IL + -tg) + U{-r - tg)\ 


e&flr 


+ ||/3||^oo I -hr 11^1100+ V{-r-ti) 

\ £eo+ e&Qr / 

« l|o||,»U. Ik'lloo + 2 ||£/|L + 2 ^ £/(nA)) + WDwJl, lIr'IL + 2 IVIL + 2 ^ l"(nA)), 


n=l 


n=l 


where we isolated the points in T that are closest to r and —r as in (79)-(80) and we 
used the fact that a regular spacing of the tg, h G hi, maximizes the sum as in (81). 
Since '^n^iU{nA) < oo and '^n^gV^nA) < oo, the Weierstrass M-test tells us that 
~ ~ converges uniformly on every compact set [—r, r], r > 0. 

Thanks to [31, Thm. V.2.14] this implies that the function A*gCo is differentiable on M, and 
that its derivative equals 


Vt G M, (^*co)'(t) = ^ {agu{t - tg) + (3gv\t - tg)) 

£&n 

for a,/3 E We next show that there exist a,/3 E such that {A*gCQ){tg) = eg, for 

all I eVL, and |(^*co)(t)| < 1 for alH G M \ T. To this end, we seek a,/3 E such that 


(^*co)'(t£) = 0, 


(93) 


for all i E Q. In developing an approach to solving the equation system (93), it will turn out 
convenient to dehne the operators 


Up: h~(fi) 

a = {oigjg^fi 


and 


V„: r 


^p. ^ ^(fl) 




-{Q) 

Y. a^u^P\tg-tY 

mGn 


E l^mV^^\tg-tm) 

m(gei 




ten 


where p E {0,1}. We defer the proof of Up and Vp, p E {0,1}, mapping £°°(r2) into £°°(r2) to 
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later. The equation system (93) can now be expressed as 


WqQ; + Vo/5 — £ 
Uia + Vi/5 = 0. 


( 94 ) 


If both Vi and Mo - VoVf'Mi are invertible, then, as in [3], one can choose a = {IAq — 
VqVi^Ui)~^£ and f3 = —Vf^Wia to satisfy (94). The Neumann expansion theorem [32, 
Thm. 1.3, p. 5] now says that ||X — (n'(0))“^Vi|| < 1 and ||X — (Wq — VoVh^Wi)|| < 1 are 
sufficient conditions for Vi and Wq —VoVf to be invertible. We next verify these conditions. 

A.2.1 Vi is invertible 

Fix a sequence /5 e £“(12), dehne ( = {I — (n'(0))“^Vi)/5, and let £ G 12. We then have 

Ce = I3i- (^''(0))■^ ^ I3mv'{ti - tm) 

= -(n'(0))^^ ^ PmV'iti - tm) 

= (oR"(0))-l l3mV\t,-tm), 

mSO\{£} 

where we used n'(0) = -i?"(0) > 0. With (92) and i?"(0) = we obtain 

2^2 
TT 


^ V —\^m\W{k-t 

< ^ TT 

2^2 
TT 


^- \\/^\\e°- Y] y{te-tm) 

TT ' ^ 

ll/3|l,» E ^ ^ 


mGQ\{i} 


T^Jc\k-tm\ 4 :a^Ic J 


(95) 

(96) 


We further upper-bound (96) using the same line of reasoning that led to (81). Specifically, 
we make use of the fact that V is non-increasing on (0, cx)) and that the minimum distance 
between points in T is A. This implies that ~ maximized for 


{t£ — tm'- m G Q \ {£}} C |nA: n G X \ {0}}. 


With (92) this gives 




( 97 ) 
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As n ^ 1, we have R{nA) = exp j ^ exp^—j, which when used ir 

to a further upper bound in terms of the following power series 


Va;e(-l,l), ln(l-a;) =, 


n=l 


X 


oo 

= nx' 


X 


(1 — xR ^ ’ 1 — X 

^ ' n=l 


all evaluated at a; = exp(—< 1. Putting things together, we obtain 


X-(n'(0))-iVi =sup 


/3^0 


11/911 




tt/cA 


In ( 1 — exp 


ttA" 


Aexp(-|J? 


+ 


4(j^ 

2exp(-^ 

TT^ 


Dehning the functions 


2^VPl-exp(-|Jf))= l-exp(-lfx) 


Vx > 0, ^{x) := — In (l — exp(—TTo;^)) 

x^ exp(—TTo;^) 


'ix > 0, il){x) : = 

Vx > 0, ^{x) = 


(1 — exp(—TTX^))^ 
exp(—TTo;^) 

1 — exp(—TTo;^) ’ 


we can then write 


X-(n'( 0 ))-'Vi ^ 


/A 




The functions and are non-increasing on (0, cx)), as their derivatives satisfy 

w n \ -27ra;exp(-7ra;^) 

^->0. l-exp(-xx^) 

Vx>0, ('(x) = _ 2xxexp(-2xx^)^ ^ ^ 


As for -ip, we hrst write 

\/x > 0, 'ijj{x) = 
and then show that the function 


1 - exp(-7ra;2) (i _ exp(-7ra;2))^ 
a;exp(—7ra;^/2)^ ^ 


1 — exp(—TTX^) / \2sinh(7ra;2/2) J 


X 


'ix > 0, rj{x) := 


X 


2sinh(7ra;^/2) 


(97) leads 

OO 

n=l 
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is non-increasing on (0, cx)) by computing its first derivative: 


\/x > 0, ri'{x) 


2sinh(7ra;^/2) — 27rx‘^ cosh(7ra;^/2) 
(2sinh(7ra;2/2))^ 

cosh(7ra;^/2) (tanh(7ra;^/2) — nx"^) ^ ^ 
2 (sinh(7ra;2/2))^ 


where the inequality is thanks to tanh(7ra;^/2) ^ 7ra;^/2 ^ ttx"^, for all a; > 0. Therefore, the 
function ip is also non-increasing on (0, cx)). Since by assumption A > 4a and A > l//c, we 
get 


X-(n'(0))-iVi ^ 


f^(2) + 2^(2) 

/cA 


+ 2^(2) 


< + 2ip{2) + 2^(2) ^ 3.71 ■ 10"^ < 1. 

71 

It therefore follows that Vi is invertible. Furthermore, according to the Neumann expansion 
theorem, the operator norm of satishes 


vr 


-II 




|t;'(0)| 


-1 




2(7^ 


||X-(r;'(0))-iVi|| ^ TT-2(^(2)+7r^(2)+7re(2))' 


(98) 


A.2.2 Uq — VqVi ^Ui is invertible 

We start by noting that thanks to the triangle inequality. 


x-(Wo-VoVriWi) < 


< 


i-Moii + ilv„ii||vr‘l!il«.il 

2^-"«ll+l_||I_(„,(0))-.V. 


(99) 


An upper bound on ||X — Uo\\ can easily be derived using arguments similar to those employed 
in Section A.2.1 to get an upper bound on ||X — (n'(0))“Wi||. Specifically, for a G £“(12), 
the sequence C, = {X — 7/o)a obeys 


E |cim| tm )l ^ ll«llo» 


mSO\{q 


^ 2||a| 


E 


for all £ G 12, where we used 


exp 


r2iGri\{£} 

7inA^\ II CK 


27r/c 


1 tra) ^ 2 llcrll^oo ^ ^ 


R{nA) 
^ 27r/cnA 


40-2 J tt/cA 


In ( 1 — exp 


ttA^ 

4(t2 


Vt G M\{0}, |M(t)| ^ 


m) 

27r/c|tr 
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This implies that 




( 100 ) 


tt/cA tt/cA tt 

Next, we compute an upper bound on ||Wi||. To this end, we £x a G £°^(Q) and set ( = Uia. 
Since -R'(O) = 0 and sinc'(O) = 0, we have n'(0) = 0, which, combined with (91) gives, for all 

£ e n. 


^ ^ ^ |cim| 1"^ (tl tm)\ ^ ^ |cim| 1"^ (ti tm)\ 

1 


^ IIq;|| 


E 

mSO\{£} 


1 1 

+ V-^ + 


4o-2/c \ti - tm\ 27r/c \te - U, 


R{tl tm) 


^ 2||a| 


^ ll«l 


= a 


4/cO-' 


^ f2 R(nA) + f2 > + f2 


n=l 


1 / nnA^\ (2 

2^exp(-^ 1 + I - + 


2/cO- 


nA ^ 2Tifcn^A? 

=1 n=l •’ 

1 


n=l 


4(7^ 


1 f TinA? 

A ■ ./.Av£n“P 


2/,a2(l-exp(-^)) VA tt/.A^ 


- ^ + 


In 1 — exp 


4(7^ 

ttA^ 


( 101 ) 


4(7^ 


where in (101) we used the fact that for 1, 1/n^ ^ 1/n and R{nA) = exp^- 
exp^—j. Based on the npper bound (101) we can now conclude that 


2^2 


40-2 






2/cO-^ ' ''\2a‘^) \A tt/cA^J ‘ 


Setting 


\/x > 0, p{x) := x^^{x) = 


x"^ exp(—Tur^) 

1 — exp(—TTo;^) ’ 


( 102 ) 


we can rewrite (102) as 




/,A2 "V2^V VA TifcA\ 

We can verify that p is non-increasing on (0, cx)), which finally yields 

2p{2) + (2 + l/7r)<^(2) 


II«.II«¥# + -/’(2)(x+ ^ 


/.A" 


A 7r/cA2 




4(7 


(103) 
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It remains to upper-bound ||Vo||. To this end, we fix G £°^(Q) and define ( = Vo/3. As 
u(0) = 0 and 

m 


VteM\{0}, |u(t)| ^ 


4a2/c’ 


we get 


/A 

EiA 

111 l^(^^ 

tm) 1 

E i'3".i 

\v{tl - tm)\ 


mGQ 





1 /3| £00 

E 

mSO\{£} 

to 1 

7^ ^ 

1 ^ 1 /^l 

2(t2/c 

OQ 

Efl(nA) 

n=l 



00 

/ ttuA^ 

A _ 

. exp(-^) 

2fca^ 

4a2 

n=l ' 

) 2/cO-2 



for all £ G This yields 


IIVoll ^ hi) = dil g M2) 


(104) 


2/cu2 2a " 2a ’ 

where the last inequality follows from A > l//c, A > 4a, and the fact that p is non-increasing 
on (0,cx)). Finally, using (98), (100), (103), and (104) in (99), we obtain 

III - W|| < + p(2) 1,(2) + (1 + l/(2.))^(2)] . 

" "" IT 2ir - 4 (s 3(2) + ir<P(2) + ir^(2)) ' ’ 

where W '.= IAq — VqVi^IAi. Again, applying the Neumann expansion theorem, we can 
conclude that the operator W =Uq — VoVf^Wi is invertible and that its inverse satisfies 

^( 2 ) p(2)[p(2) + (l + l/(27r))v^(2)]V' 


W 


-II 




1-llX-Wll 


^ 1 - 


71 271 - 4{(p{2) + 7r'ijj{2) + 71^(2)) J 


(105) 


For later use, we record that for the choices a = {Uq — VqVi ^Ui) and [3 = —Vi ^Uia, we 
have 


«||,.o ^ llw^ll I|£||,^ = llw^ll 

< A _ ^ _ P(2) W) + (1 + l/(27r))^(2)] 
77 27r - 4 ( 99 ( 2 )-h 7r'?/>(2)-h 7r,^(2)) 


-1 

< 1.01 


(106) 


and 




K^\\\m ii« 



< 

’ ^ TT - 2(99(2) -h 7r^/>(2) -h 7r^(2)) 
p(2) [p(2) + (l + l/(27r))^(2)] \-^ 
27r - 4(99(2) -h 7rilj{2) + 77^(2)) J 


2p(2) + (2 + l/7r)(p(2) 

4a 

^ 5.73 ■ 10“V. (107) 
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In the remainder of the proof, we exclusively consider cq with a = (Wq — VoV^ ^Ui) and 

A.3 |(A*co)(^)| < 1 for all t G M\r 


A.3.1 


(Gl*co)(t)| < 1 for all t G M\lj£6r 



1 

7/c 


, ti + 


1 

7/c 


Take £0 G and let £1 G hi be the index of the point in T that is closest to and satishes 
ill > iio- Take t G and note that the interval is 

non-empty because — tii\ ^ A > T > 7T- Without loss of generality, we assume that 
1^ — iiol ^ \t ~ iiil^ which implies \t — till ^ \tii — t^l /2 ^ A/2. We set h := \t — and 
note that h ^ The following holds 

' Jc 


||^||_ 2^ + 11/911 _ 2^ + ||,,||_ + II/9II 

TrfrA 


left 


-2k).iii«ii ihhiiuii ^(t), I..,I 

2iTfJi ^ 4cVo " ir/„A ^ 

f \\n,\\ R{i-ii) II oil R{i - i- 

' ^ 27r/e |f - ^ 4a2/, 

teo\{£o/i} ^ 

„ „ i?(h) i?(h) „ „ i?(|) 

n ^(77A) i?(nA)\ 

+ 2 g (^||a||,oo + W^Wi^ j 

r(w) r(-) 

11 


tt/cA / 

+ 11/^11 


ll/^ll^oo / 

fc / 

exp(-i) . 

a 


+ ll/^llf°o 


exp(-7r; , 

~r ^ poo 

a TT 


^ + Wh¬ 


im 

(109) 

( 110 ) 
( 111 ) 


where (108) and (109) follow from h ^ and (110) and (111) can be derived invoking the 
assumptions A > l//c and cr = W_ 
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A.4 A*co is concave on te - 


7fc 


Let £ G f2. We show that t h-)■ A{t) := |(A*co)(t)| is strictly concave on 
Since |i?|, |sinc|, and |sinc'| are all symmetric, A is symmetric as wel 


ti , te + 7^.^ 

, and therefore, it 


suffices to show that A''(t) < 0 for t G te,te + ^ ■ Since \e(\ = 1, we can write 


Vt G M, A{t)= {A*Co){t) = 


1 _ 

(A*Co)(t) 


££ 


where A^it) := Re 
A{t) ^ 0} we have 


{A*co){t) 




, Ai{t) := Im 


{A*gC0){t) \ 


££ 


= V{MtW + {A,{tW 

|, for t G M. With A := {t G M: 


Vt G A, A'\t) = 


AUt)Anit) + A"j{t)Ai{t) + l(Al*co)'(t)l {A'^{t)An{t) + A'j{t)Ai{t)f 


For A to be concave on 

'it G 


A{t) 

te,te + , it therefore suffices to show that 






) + A/(t)A/(t) + (A*gCoy{t) < 0. 


Let t G 


ti, te+ ■ We have the following 


— }v(t- tm) 


^R^t) = ^ ^Re I ^ I - tm) + Re 

lit — t()+ 'Re{ — \ v{t — t() + I^Re I I 


= 'R.e { — )■ u{ 
£e 




>Rer^ 


+ Re 1^1 n(t - tm) 

- ||/3||^oo |n(t-ff)| - ^ (^\\a\\f,^\u{t-tm)\ 

A ™^o\ ( 

+ WW^oo \v{t - tm)\y 

With a = = e — {I — W“^)e, it follows that 


mGft\{i} 


Id 


Re = 1 - Re 


[(X - 

£e 


^ 1 - 


[(X - W-^)e], 


£e 


^ 1 - ||X-w^ll 

= 1-1|>V^(X->V)|| ^ 1 - IIX-WII ^ 0.999998, 


(112) 
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where (112) is due to (105). Next, it follows from t — te ^ and c = ^ that 


4/c 


u{t - te) = R{t - ti) sinc(27r/c(t - U)) = exp( ) sinc(27r/c(t - te)) 


^ exp 


4a^ 

TT 


49A-4.V-^(t) 


dTT 


27r 


= exp — — sine, „ 

49/ V 7 

Since |sinc| ^ 1, we have 

||//||^oo iv(t -te)\ ^ ||//||^oo - t) \ |sinc(27r/c(t - te))\ ^ ||//||^oo \R'{t - te)\. 


( 113 ) 


As \R'\ has its maxima at the points ■R(T\j\ with corresponding maximum values 
4exp(-l) we get 


||//||^oo |n(t -t^)! ^ 


m, 


, 1 \ /TT 

a ""n 2/V2 


(114) 


As for every m G \ {£}, we have |t — tm\ ^ \ti — tm\ — \t — 'ti\ ^ ^ h holds that 


^ ^ ^ tyTj)! H“ 11/311^00 ^ ^ ||o; 

OO / 

+ 2 f lla 


7i?(y r[^ 

” +\m,: 


mSO\{£} 


127r 


4aVc 


n=l 




27r/cnA 


4a2/c 


^ lla 


7exp(-^(f)") 


+ 


11/91 


12 


^ exp -TTlyl +||a||,^ ^ 


127r 

a(2) , II//II, 


+ 


(j 


2^(2). 


Combining (106), (107), (112), (113), (114), and (115) yields 

AR{t) ^ 0.673. 


(115) 


(116) 
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Next, we derive an upper bound on A'^(t): 


men ^ Km) 


— \ u'{t - tm) + Re < — ^ v"{t - tm) 


/ 

\ J 


^ Re <j — !> u"(t - ti) + ||/3||^oo |n"(t - t()\ 


( ||a||^oo \u''{t - t^)\ + \\I3\\^^ \v”{t - tm)| )• 

mSO\{£} 


For all t G M, we have 


u”(t) = R”(t) sinc(27r/ct) + 47rfcR'{t) sinc'(27r/ct) + {27ifcYR(t) sine”{2nfct). (117) 


The function t h-)■ R”{t) sinc(27r/ct) is non-decreasing on [O, , since, on this interval, R” is 

negative and non-decreasing and t h-)■ sinc(27r/ct) is positive and non-increasing. The function 
t H-)■ 4:7ifcR'(t) sine'{2nfet) is non-decreasing on [O, , as both R' and t h-)■ sinc'(27r/ct) are 

negative and non-increasing on this interval. The function t h-)■ (27r/c)^R(t) sine”{2 t[ fR) 
is non-decreasing on [O, , as, on this interval, R is positive and non-increasing, and 

t H-)■ sine”{2 t[ fct) is negative and non-decreasing. Taken together, it follows that u" is non¬ 
decreasing on [O, . Since t — t( G [O, , we then have 


u”it-te) ^ u”{ = fl 


„ , Stt \ . (2tk 

Stt -1 sine — 

49 / V 7 


327r2 


, A Z ■ //I 

-|-47r sine — 
7 


47r\ 


2ir 

Sine' , ^ 

7 V 7 


where we used a = ^. Combined with (112), this yields 


Re <( — [■ u”(t — ti) ^ —6.46/^ 

££ 


Since 


,, . // X coslx) sinla;) 

\/x G ]R\{0}, sine (x) =-^— and 


sinc'Ta;) = 


sin(a;) 2cos(a;) 2sin(a;) 


X 


x^ 


X 


x^ 


X-^ 
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we get the following from (117): 


KWK I i + 


20-2 40-4 j 271 fc \t\ 


R{t) + ATrfc^^^^ 


+ 


2a2 \27rfc\t\ (27r/c)2 |f| 


m 




+ 


+ 


27r/,|t| ' (27r/,)2|t|2 (277f,)^\tf) 


m 


71 \t\ 71 


S/cCT^ 0-2 


— + 271 fc + 


1 2 
TTT + -TTo + 


4c^Vc/ |t| \tf Tlfcltf] 


R{t). 


As a result, we have the following chain of inequalities 


\u'{t-tm)\ ^ 2 y^ 

meO\{£} n=l 

TrnA^ 


nA 71 


7T 


S/cCT^ 0-2 


+ + 27r/c + 733^ - 7 - + 


+ 


nA n 2 A 2 TlfcTl^A^ 


R{nA) 


ttA ^ 

^ —r— 2^ n exp 


4o-Vc 


n=l 


4(t2 


27r ^ 

n=l 


TinA? 

4(t2 



We can now dehne 77 ( 0 ;) := x^'ip{x), for all a; > 0, which can be shown to be non-increasing 
on (0,cx)). This yields 


mGQ\{£} 


< 


7T 




2 . /A 
^ 2 ^V 2 a 


4 / A\ 2 


47r/c 

^ + 2 a 2 /,A 


167r77(2) + 327r^(2) + Utt + 


+ - ) a( 2 ) 

TT 


fl 



Combined with (106), we get 

l|a||£- \u{t -tm)\^ 1.196 ■ 10-Vc- 

77iGr2\{£} 
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We have, for all t G M, that 


v"{t) = sinc(27r/ct) — ATrfcR''{—t) sine'{2 tt fj:) + {2 tt fc)"^R'{—t) sine"{2 tt fet). 


Therefore, we get 




4a4 2a^J \27rf,\t\ (27r/,)2 |t|V 

(n f \2 ^ \ \ f ^ 2 2 

+ ^tc) ^ \t\ 1^27r/,|t| + (27r/,)2 \tf + (27r/,)3 |t|3 

/ TT^ \t\^ 7^ StT ^ 1 \ 

^ I 16a6/, + 2a4 + 8f,a^ + |t| + Itf ’ 


which leads to the following chain of inequalities 

m ^ TT^nA Stt tt^/c 27r 1 \ r-./ 

|t^ (tf-U)l ^ y iQ^Gf + 2^4 ^ 8/ca4 ^ ^ ^ a2nA ^ ^ 

meO\{£} n=l jc / 


tt^A^ “ 2 f ^nA‘^\ tt^A^ 

•’ n=l ^ ' n=l 


nexp 


( 118 ) 


OTT V "v 

n=l 

47r ^ 1 


7rnA^\ 27r^/c 


1 / 7rnA^\ 1 1 / TrnA^ 

n=l ^ \ 


Now, in (118), we recognize the power series 


Va;e(-l,l), 


2 „ x(x + 1) 
nX" = -1 -^ 

(1 — XR 


evaluated at a; = exp^—which leads us to set 


\/x > 0 , r(a;) := 


exp(—7ra;^)(exp(—TTo;^) + 1) 
(1 — exp(—7ra;2))3 
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This yields 


I /// I f f Stt 27r' 

-o\r«i J ^ \ / \ / \jt- 


mSO\{£} 


(T^ 


2 a / 


47r 1 

+ 




A 


(T^A /cU^A^y ^2(7 

TT^ f 27r^ f 

^ —^71 1 + 77yyx7 + 


2/ca4 V2a/ da^A \2a J \4:fc<J^ 


5. ,A4V(7) 


47r 1 

+ 




A\ 


- 


< 


a2A /,a2A2 J ^ \2a J 

1287r27(2) + 647r27(2) + (3207r + S27r^)^{2) + (647r + 16)(p{2)\f^, 


where we set 7 ( 0 ;) := x^V[x), for a; > 0 , and used the fact that 7 is non-increasing on ( 0 , cx)). 
Combined with (107), this results in 


WWt^ Y. 8.34 ■ 10-Vc. (119) 

Finally, we have A"^{t) ^ —22.1/2. Multiplying (119) with (116) leads to Aji{t)A'l^{t) ^ 
— 14.6//. Exactly the same line of reasoning can be applied to get Ai(t)A'l(t) ^ —14.6//, 
and therefore, 

AR{t)A!'^{t) + Ai{t)A!',{t) < -29.1//. (120) 

I 12 

It remains to find an upper bound on |(^*co)'(t)| . We have 
|(7l*Co)'(t)| ^ X] ( ll«ll£°° (^ -tm)\ + ll/^ll^oo \v'{t 

^ llttll^oo \u'{t -t()\ + ||/S||^oo \v'{t - ti>)\ (121) 

+ Y ( ll“il£°° “ ^m)| + ll/^ll^oo |n'(t - fm)| )• (122) 

mSO\{£} 

We can derive upper bounds for the terms in (121) by noting that 

l«V)l < ..'(^) = smc(7) smc' (7) ( 123 ) 

^e.p(-7).,.(-),2.7exp(-|)s.„.(7) 

and 

|F(t)|^F(0) = -i?"(0) = -^ (124) 
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for all t G 0; ^ • Indeed, we have seen that u"(t) ^ 0 for all t G 0, ^ , which implies that 
u' is non-increasing on 0, ^ . As m'(0 ) = 0 , this means that u' is non-positive on 0, . 

Therefore, \u'\ is non-decreasing on 0, ^ , which results in (123). The inequality in (124) 
follows from the fact that In'I is decreasing on 0, ^ , as we show next. We have 

|_ ' /c J 

Vt G M, v{t) = —R”{—t) sinc(27r/ct) -|- 2ti fcB!{—t) sinc'(27r/ct) 

= —R"{t) sinc(27r/ct) — 27r fcR'{t) sine {2%fct). 

As the functions t h-)■ R''(t) sinc(27r/ct) and t h-)■ 27ifcR'(t) sinc'(27r/ct) were shown to both be 
non-decreasing on 0, , we get that v' is non-increasing on 0, ^ . Moreover, we have 

'Jc ijc 

”'(t) " 

Hence, v' is non-negative on 0, This allows us to conclude that |n'| is non-increasing 

on 0, , which establishes (124). It remains to upper-bound the term in (122), which is 

done as follows: 


( llall^oo \u{t -tm)\ + \v'{t - tm)| ) 

'fTiGn\{£} 

^ ||a||foo 17(t - tm)||/3||^oo H(t - tm)) 

„( \ . .yet, r2p(2) + (2 + l/!rV(2) 

^ 11 Cl 11 £oo I l+||/3||^ooHI^^ 1 T||ci||^oo ^ 


^ 0 ^( 2 )+ 2 ^( 2 )+ 2 ^( 2 )) 


< ll«| 


1 7/ 49/, 

4a2/c ^ 6 ^ 727r 




7 TT Stt 
^ ^ ^ 2 W? 




+ i|ci||^oo [2p(2) -I- (2 -I- l/7r)v9(2)_ fc H- j—(^(p{2) + T[ip{2) + 7r.^(2)) 


«IMI,«(4+^ + ^)/.exp(- 


5767r\ , \\I3\\,^ /7 


X + 27r+ — /cexp - 


A 4.05/e. 


+ ||a||,» [2p(2) + (2 + l/7r)(p(2)]/e+ ^^(4(p(2) +4x^(2) +47r/(2))/e 

(125) 
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Putting (120) and (125) together yields 

AR{t)A"^{t) + Aj{t)A",{t) + |(^*Co)'(t)f ^ -12.68/2 < 0, 
which completes the proof. 

B Proof of Theorem 11 

We could prove Theorem 11 following similar arguments as in the proof of Theorem 10, 
namely by choosing a function cq G L°°(T x Z) of the form 

L 

Vr e T, VA; e {-iP,..., iP}, co(r, k) = Y, - r)e-2"'=*^ + p,g'{te - 

e=i 

and determining a := and f3 := such that the uniqueness conditions (41) and 

(42) are met. It turns out, however, that a more direct path is possible, namely by choosing 
a function cq G L°°(T x Z) of slightly different form and then reducing to a case already 
treated in the proof of Theorem 10; this approach leads to a substantially shorter proof. We 
start by dehning this function cq G L°°(T x Z) as 

Vt e T, Vfc e {-A-,..., K}, c„(t, k) := ^ (a,p(r-i,)e-““<+ftg(T-«,)e-='““), 

where p: T —)■ C and g: T —)■ C are dehned (for reasons that will become clear later) as 
Pir) := ^ 

nGZ nGZ 

for r G T, with 

_ ri/a 

Pn ■= y/2aexp (^—27ra'^n'^) / exp (—exp (—STra^nn) du 

J-l/2 

ri/2 

g„ := —27ri(TV^^exp iy—2Tia‘^n^') / (u + n) exp (—exp (—STra^nn) dn, 

/-1/2 

for n G Z. We hrst verify that the resulting function cq is, indeed, in L°°(T x Z). This is 
accomplished by showing that the functions p and q are well-dehned and are in L°°(T), that 
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is, by verifying that b«l < Xlnez Indeed, we have 

j‘l/2 

\pn\ = v^^y~~^exp (—27r(T^n^) / exp (—dTrcr^M^) exp (^—Sna'^nuj du 
nez nez J-1/2 

^ j exp (—27r(T^(4M + n)n) dw (126) 

nez -^-1/2 

^ /-i/a 

= C + / ) / exp (—27r(T^(4M + n)n) dw, (127) 

nez -^-1/2 
|n|^3 

where (126) follows from exp(—dTra^ti^) ^ 1, for all u G [—1/2,1/2], and we set 

2^ /.1/2 

C := / exp (—27r(T^(4M + n)n) dw < cx). 

n=-2'^“l/2 

To see that the sum in (127) is hnite, hrst note that for n ^ 3 and u G [—1/2,1/2], we have 

(4m + n)n = (4m + n)|n| ^ (—2 + n)|n| ^ \n\. 

Similarly, for n ^ —3, we get 

(4m + n)n = —(4m + n)|n| ^ —(2 + n)|n| ^ |n| . 


It therefore follows that 


7^/2 rl /2 

yj / exp (— 27 r(j^( 4 M + n)n) dM ^ / exp(— 27 r(j^|n|)dr 

J-l/2 „^.y, J-1/2 


oo 

= 2 exp(—27r(7^n) < cx). 

n=3 
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This concludes the proof of \Pn\ < oo- Similar reasoning shows that \^n\ < oo. 

For t G T, we then have 


- 1/2 


i^*gCo){t) = ^ [ co{T,k)g{t-T)e™dT 


(128) 


2F:+ 1 


t=i 


r>l/2 

I / Pij - - 'r)dr ) DK{t - te) 

- 1/2 

n /2 


1 


2f: + 1 


+f^e \^J q{r - te)g{t - r)dr ) Dxit - ti) 

L 

(^oi£P{t — ti)DK{t — t^) + PiQ{t — te)DK{t — tg)^, 

e=i 


where Dk is the Dirichlet kernel, that is, 


Vt G T, — 


sin((2iF + l)7rt) 
sin(7rt) 


and P and Q designate the cross-correlation between the functions p and g, and q and g, 
respectively, that is. 


Vt G T, 
Vf G T, 



P{r)g{t 

(l{r)g{t 


T)dr 


r)dr. 


Note that since g and r h-)■ co(r, k), k G {—K,... ,iF}, are all 1-periodic, we can integrate 
over the interval [—1/2,1/2] in (128) (instead of [0,1] as done in (18)) and in the remainder 
of the proof. We next derive an alternative expression for the function P. As in (54), we 
have 

VfGT, (/(t) = 5^<?ne“, 

nGZ 

where gn '■= v^^exp(—27r(T^n^), for all n G h. The nth Fourier series coefficient of P is 
then given by Pnfi'n, and we show that the Fourier series Xlnezconverges to P(t) 
for all t G T using Dirichlet’s theorem [33, Thm. 2.1], whose applicability conditions we 
verify next. Since Yhn&z \Pn\ < Snez ^ |g 2 ™t| _ f G T, by 

the Weierstrass M-test, the series and Xlnez converge absolutely and 

uniformly. This implies that the functions p and g are both continuous on T. Moreover, 
g is continuously differentiable on M as \^9n\ < oo. As a result, the function P is 


61 



continuously differentiable on M, and by application of Dirichlet’s theorem, it follows that 


Vt e T, P{t) = ^ QnPn 


For n G Z, we have 

ri/2 

gnPn = V^exp (—27r(T^n^) v^^exp (—27r(T^n^) / exp (—exp (—STra^nw) du 

J - 1/2 

H/2 

= 2(Texp (—dTra^n^) / exp (—exp (—STra^nti) dw 

J - 1/2 

rl/2 

= 2a exp (—47r(j^(M + n)^) dw. 

J-l/2 

Now £x t G [0,1). If t = 0, we have 

l‘l/2 

P{t) = -P(O) = gnPn = 2a / exp (—47ra^(M + n)^) dw 

neZ 

/ n+1/2 1*00 

exp(—47ra^n^)dn = 2a exp(—47ra^n^)dn = 1, 

-1/2 J-oo 


and if t ^ 0 , we get 


P{t) = ( 20 '/ exp (— 47 ra^(M + n)^) dw | e 

nez V -^- 1/2 / 

/‘I /2 ^ 

= / 2 a V exp (- 47 ra 2 (M + n)^) e“dM 
/‘I /2 

= / 2 a V exp (- 47 ra 2 (M + n)^) 

— 1/2 .^^1? 


I 2a V exp (- 47 ra 2 (M + nf) 
“V 2 
/‘ 1/2 

1 //'(m)(P(—M) dM = .^(0). 

- 1/2 


(129) 

(130) 

(131) 


Here, ip is the 1-periodic function defined by ip{u) := u G [—1/2,1/2), and tjj and ^ are 

given by 


Vm G T, 

il){u) 

:= 2 a exp (- 47 ra 2 (M + nf) 

nGZ 

1 . 1/2 

'ix G T, 

i{x) 

:= / (p{u)ijj{x — u)du. 

J-1/2 
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The order of summation and integration in (129) is interchangeable thanks to 


/ i/^ _ rn+i/z 

exp (—47r(T^(M + n)^) du = / exp (—du 

-1/2 in-1/2 

/ oo 

exp (—du < cx). 


The function ^ can be expanded into a Fourier series. Specihcally, it holds that 

VxeT, e(a:)=5^</^ni’ne2™", 


neZ 


where and denote the nth Fourier series coefficients of (p and //>, respectively. We have 


pll2 

cp^= / (^(a;)e-2™^da; = 

i-l/2 


(-l)"sin(7rf) 


/■ 1/2 

/ = 

'-1/2 7r(t - n) 


and 


iJn= [ i/’(a;)e-2™"da; = [ 2u ^ exp (-47ra2(a; + m)^) 

>'-1/2 i-l/2 

^m+1/2 

= 2a^ / exp (-47ra2u2) 
mez -'"^-1/2 


= 2a exp(—47r(T^u^)e^’^*‘'* 

J —OO 

7r(t — n)^' 


= exp 


4(t2 


for n G Z. It follows that 


VteT, = = 


nGZ 




Ti{t — n) 


We then get 






Ti{t — n) 


n(t — n)‘ 
4(7^ 


n(t — nY 
4(7^ 


{2K + 1) ^^sinc((2iF + l)7r(t — n))R{t — n), 


nGZ 


where R was dehned in (15). Similarly, we can show that 

Vf e T, Q{t)DK{t) = {2K + 1) 5]] sinc((2iF + l)7r(t - n))R!{n - t). 
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This finally yields 


L 

VteT, {A*gCo){t) = ai sinc((2if + l)7r(t — te — n))R{t — te — n) 

£=1 nGZ 

+ /Si sinc((2i^ + l)7r(t — — n))R'{n — t + tg)^ 

L 

aguit — tg — n) + (Sgv{t — tg — n) j, 

£=1 n<^'L 

where we set 


'it G M, u{t) := -R(t) sinc(27r/^t) (132) 

Tit 

it G M, v(t) := -R'(—t) sinc(27r/'t) = —sinc(27r/^t) (133) 


as in (86) and (87) with f'^ := K + 1/2. Analogously to the proof of Theorem 10 we can 
define the operators 


Up-. 




a = {agjf^^ 


X) amU^P\tg - tr 

m=l ngZ 


ni 


e=i 


and 

Vp: 

= {/^ejeen 

where p G {0,1}. Then, given e 
equation system 


E - t^ - u) 


m=l 


£=l 


with le^l = 1, £ G {1, 2,..., L}, we can solve the 


UqCX + Vo/3 — E 

Uia + Vi/3 = 0 


(134) 


to determine a G and /3 G such that the interpolation conditions {A*gCo)(tg) = Eg, for 
all £ G {1,2,..., L}, are satisfied and Vl*Co has a local extremum at every tg, £ G (1, 2,..., L}. 
As in the proof of Theorem 10, if the operators Vi and Uq — VoVf are invertible, then 
one can choose a = {Uq — VoVf^7/i)“^e and /S = —V^^Uga to satisfy (134). Proving the 
invertibility of Vi and Uq — VoVf is essentially identical to the corresponding part in the 
proof of Theorem 10 with fc replaced by f£ Verifying that |(^*Co)(t)| < 1 for alH G T \ T, 
where T = {tg}f^i, is also done in a fashion similar to the proof of Theorem 10 (see Section 
A.3). 
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